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Summary 


The seasonal variation of the pack ice surrounding the Antarctic 
continent has a great effect on the levels of microseismic activity at Scott 
Base. Microseismic storms are in general of two types, one with periods 
1-3°5s, the other with periods 4-10s. ‘The short-period activity is due 
to events within the Ross Sea and is at a maximum during January and 
February when the Ross Sea is clear of ice. The long-period activity has 
a maximum during March and April and this is thought to be due to 
swell from peripheral cyclones penetrating the Ross Sea and creating 
microseisms after passing the continental shelf. 

It is found that the microseisms of any period (7) have a limiting 
amplitude (A max) and that for the short-period microseisms Amax « T?” 
while for the long-period microseisms Amax « 76°. The difference be- 
tween these two relationships can be satisfactorily explained on the 
basis of frequency selective decay of the waves causing the long-period 
microseisms as these waves travel from the cyclone area to the continen- 
tal shelf. ‘The fourth-power relationship between maximum amplitude 
and period is similar to that obtained for single storms in America by 
Romney. It is suggested that this relationship might provide a useful 
test for a quantitative theory of microseisms, such as the standing wave 
theory, or alternatively that it can be used to investigate transmission 
and absorption processes in the crust. 

The data appear to confirm the requirement of sea roughness for the 
generation of microseisms, the half-period relationship between micro- 
seism and sea wave periods, and the role of the continental shelf and 
other crustal features as barriers to microseisms of periods less than 
738. 
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microseisms as these waves travel from the cyclone area to the continen- 
tal shelf. The fourth-power relationship between maximum amplitude 
and period is similar to that obtained for single storms in America by 
Romney. It is suggested that this relationship might provide a useful 
test for a quantitative theory of microseisms, such as the standing wave 
theory, or alternatively that it can be used to investigate transmission 
and absorption processes in the crust. 

The data appear to confirm the requirement of sea roughness for the 
generation of microseisms, the half-period relationship between micro- 
seism and sea wave periods, and the role of the continental shelf and 
other crustal features as barriers to microseisms of periods less than 
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Introduction 


In 1957 March a seismological station was established at Scott Base (77° 51'S, 
166° 48’ E) as part of the New Zealand I.G.Y. Antarctic programme. The seismo- 
meters were three single component Benioff reluctance-type instruments of 
nominal 1 second period. The short-period coils fed three 0-2s galvanometers; 
the long-period coils on the vertical and east-west instruments were connected to 
25s galvanometers and on the north-south seismometer to a 10s galvanometer. 
The system was not designed for the optimum recording of microseisms. Never- 
theless due to the geographical position of Scott Base and the effect of the Ant- 
arctic pack ice and sea ice certain simple well defined patterns emerge concerning 
the behaviour of microseisms. Deacon (1953) has called for the study of micro- 
seisms where “‘the phenomenon . . . appeared to be simplest’; the control of the 
sea surface effected on the one hand by the ice of the Ross Sea, and on the other by 
the Antarctic pack ice has provided some extraordinarily simplified physical con- 
ditions which do not obtain elsewhere. Reports of microseisms at stations on the 
Antarctic continent have been made before (Thompson 1950, Imbert 1954) but 
for reasons which will become obvious in the main text Scott Base has certain 
geographical advantages over these other stations. 

Microseisms of all periods from } to 10s were observed at Scott Base but in 
general they belong to two well-defined groups. The first group consists of micro- 
seisms of periods between 1 and 3:5s, and will be referred to as the short- period 
microseisms. ‘The second group has periods between 4 and 10s and these will be 
called the long-period microseisms. ‘The morphology of each of these groups will 
be separately described in this paper, and the results discussed jointly to determine 
the contributions that the measurements make to an increased understanding of 
the generation and transmission of microseisms. 


2. The location of Scott Base 


The seismological station at Scott Base lies on frozen basaltic debris at Pram 
Point on the S.E. side of the Hut Point Peninsula on Ross Island. Ross Island 
lies on the S.W. corner of the Ross Sea, and the relation of Scott Base to the 
continent, the 1000m ocean depth contour (taken as the continental shelf), the 
oceans and the Ross Sea is shown in Figure 1. This figure also shows the distribu- 
tion of ice in the vicinity of the Ross Sea at various times of the year. The maps of 
the ice distribution are drawn from material taken from a number of sources. 
These include reports of early expeditions, the U.S. Navy Hydrographic Office 
Technical Reports on Operation Deepfreeze, the U.S. Navy sailing directions for 
Antarctica and a survey by Herdman (1953). 


3. Measurement of the microseismic parameters 


In this study only the period and amplitude of the microseisms are considered. 
‘These parameters have been measured according to the procedures laid down for 
the reduction of I.G.Y. microseismic data (Rothé 1957). The raw amplitude has 
been derived by taking the mean value of the maximum amplitudes of five groups 
of the most important waves occurring in a time interval of 20min, the interval 
being symmetrical about the hours 0, 6, 12 and 18 G.M.T. The period has been 
derived by taking the mean values of the periods of the five groups of waves used 
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for the derivation of the amplitudes. Where five groups have not been available 
a lesser number has occasionally been used. 

During the past few years, following the work of Barber & Ursell (1948) 
eceanographers have discarded single-valued periods and amplitudes in favour 
of a description of the sea in terms of its Fourier amplitude spectrum. Such 
techniques are now applied to microseismic investigations (Darbyshire 1954). In 


(a) August (b) December 



































(c) Late February. (d) Late April. 
Fic. 1.—Distribution of sea ice over the Ross Sea and S.W. Pacific Ocean 
at various seasons of the year. 
— ——~-— Continental Shelf. 


the present study it has been impossible to analyse 18 months’ records in the 
manner required to give a complete description of the microseismic spectra. ‘The 
spectra of the short-period microseisms vary from the fairly narrow (Figure 2a) 
to the very broad (Figure 2b). Those of the long-period microseisms are usually 
very narrow (Figure 3). However, the periods and amplitudes of the highest 
amplitude waves are identifiable and these are considered to be suitable parameters 
for a discussion of the seasonal behaviour of microseisms. Figure 4 shows the 
period response curves for the short-period vertical component from which the 
short-period microseisms were read, and the north-south long-period component 
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from which the long-period microseisms were read. Romney (1953) has pointed 
out that the short-period Benioff seismograph when critically damped has a res- 
ponse which decreases nearly as the cube of the period for periods longer than 
4s. This is a useful attribute for Romney found, and it will again be demonstrated 
here, that the amplitudes of microseisms increase as about the fourth power of 
the period. Consequently all ground components in the 1-4s frequency range 
should be recorded equally well. 
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Fic. 2.—Short-period microseisms showing (a) typical narrow spectrum 


storm, 1958 February 24 and (b) typical wide spectrum storm 
1958 January 29. 


4. The seasonal variation in amplitude 


The short-period microseisms. The weekly means of the six-hourly amplitudes 
of the short-period microseisms are shown in Figure 5(a) for the whole of I1.G.Y. 
During both southern winters microseismic activity was virtually absent at the 
sensitivities used. (That is, ground motion of 1s period was less than 0-003 um). 
In both 1957 and 1958 the activity, which was sometimes present at low amplitudes 
in September and October, began to increase sharply during late November. The 
activity reached a maximum between late January and early March, but after the 
latter date decay was rapid and by late 1958 April the short-period activity had 
almost settled down to winter level. 

This pronounced seasonal variation must be related to the ice conditions on 
the Ross Sea. The ice begins to form in the southern part of the sea and around 
the coastline during mid-March, and by mid-winter freezing of the Ross Sea is 
probably complete though, due to storm and tide action, this ice cover may be 
broken. The ice cover extends well to the north of the continental shelf (Figure 1) 
and surface motion of the sea over the whole continental shelf is greatly inhibited. 
During spring the ice begins to break up and by mid December there is little ice 
left on the Ross Sea although the entrance to the sea is blocked to a much later date 
by broken pack ice. Thus the Ross Sea is isolated in spring and early summer. 
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The depth of water is about 500m and fairly uniform, while the dimensions of the 
sea are roughly 700km x sookm. A useful closed system appears to exist in which 
to investigate the interactions of weather, sea and crust which produce microseisms ; 
generation of the microseisms takes place directly in the continental crust which 
again provides a simplifying factor rarely encountered. 
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. 3.—(a) Typical long-period microseisms, 1958 April 20, (b) spec- 
trum of zero line crossover times, 1512-1514, 1958 April 20. 


For several years now most workers in the field of microseisms have accepted 
the concept that microseismic energy is derived from sea swell. The absence 
of microseismic activity when the Ross Sea is covered with ice is a suitable demon- 
stration that direct transmission of air pressure fluctuations to the sea bottom does 
not cause microseisms, for such energy should still be transmitted if the sea carries 
a solid floating body of ice. However since the effect of sea ice is merely to inhibit 
surface irregularities on the sea there still remains the possibility that both sea 
roughness and air movement are required to generate microseisms. 
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The long-period microseisms. Figure 5(b) shows the daily amplitudes of the 
long-period microseisms expressed as the average of the four six-hourly measure- 
ments of amplitude. The daily mean period of the most prominent activity is 
shown in the upper part of this figure. In general the long-period microseisms 
appeared to be comparatively monotonic. The exceptions are shown by a linc 
covering the period band observed or, where even this is difficult to determine, 
by the letter M. Q in the amplitude section of Figure 5(b) represents quiet con- 
ditions where microseismic activity is too weak to measure. S indicates that short- 
period storm activity obscures the long-period microseisms. 
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Fic. 4.—Period response curves of (a) the short-period vertical compo- 
nent, and (b) the long-period N-S component. 


As in the case of the short-period microseisms there is a strong seasonal varia- 
tion. Increase of activity above the winter level occurs during January and there- 
after successive storms are of increasingly greater amplitude until mid-March. 
From this time until late April the activity is at its highest level, and then the ampli- 
tudes again decrease, reaching the winter level by mid-June. Thus there is a 
difference of two months between the maximum activity of the short-period 
microseisms (January-March) and the maximum activity of the long-period 
microseisnis (March-May). This is explained in terms of the behaviour of the 
pack ice at the entrance to the Ross Sea. In Figure 1 it can be seen that this pack 
ice persists until February even though the Ross Sea is clear. By mid-March it 
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has disappeared and does not build up again until pronounced freezing takes place, 
or until the sea ice formed to the south is broken up and carried northwards by 
the winter storms. 
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Fic. 5.—Seasonal variation in amplitude of (a) the short-period micro- 
seisms, and (b) the long-period microseisms. Mean daily periods are 
shown for the long-period microseisms. Q = quiet microseismic con- 
ditions. S = short-period activity masking long-period microseisms. 


During summer it is found that the easiest passage into the Ross Sea is obtained 
by moving through the pack at about 180° longitude. Figure 6 shows the northern 
and southern limits of 47 transits of the pack ice made between 172° E and 172° W. 
This figure illustrates well the conditions which allow the generation of micro- 
seisms. 
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Fic. 6.—Northern and southern limits of the antarctic pack-ice observed 
during transits made between 172°E and 172° W. 


The early development of open sea north of Scott Base is illustrated by the retreat 
northward of the southern edge of the ice. Extensive open water exists by Novem- 
ber and this coincides with the appearance of the short-period microseisms. On 
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the other hand the northern edge of the pack ice does not move southwards to any 
extent until January. By late February little pack ice is left and this is in the region 
of the continental shelf. Refreezing of the Ross Sea commences from the shore 
about mid-March. ‘Thus the very high microseismic amplitudes in the period 
range 4~-10s are associated with the absence of ice in the northern part of the Ross 
Sea during March, April and early May. The presence of ice in this area can (a) 
inhibit the production of sea waves by cyclones or, (b) attenuate greatly the swell 
approaching the Ross Sea from distant storms. It is certainly true that much of 
the water under typical cyclones to the north of the Ross Sea can be covered by 
ice. Under winter conditions almost all of the southern half of a cyclone such as 
that shown in Figure 16 would be over ice-covered water. However alternative 
(b) involving the role of the continental shelf has been preferred for the following 
reasons. Firstly, large microseismic storms do not appear on Scott and Hallett 
seismograms until the autumn, though by the end of December the ice-edge is 
usually at 67°S latitude, which allows seas under the cylones to develop in a normal 
manner. Secondly, the amplitudes of the long-period microseisms decrease rapidly 
in late April and early May. Since the short-period activity at Hallett Station at 
the northern edge of the Ross Sea ended about April 16 in 1958, it may be presumed 
that the freezing of the Ross Sea was complete and that swell entering the Ross 
Sea after this date was subject to the damping effect of the ice. It would be con- 
siderably later, probably July, before part of the sea under the cyclones was covered 
by ice. 


Long-period microseisms of relatively small amplitudes are observed on 
occasions throughout the winter at Scott Base. It can be seen from Figure 5 that 
the periods of these winter storms are higher than the periods of the summer 


storms. 


5- The period-amplitude relationship 


The short-period microseisms.—In Figure 7 are shown the individual values of 
amplitude and period of the dominant short-period microseisms for each six- 
hourly observation during January, February and early March 1958, when it is 
assumed that the fetch was approximately constant and at a maximum. It is 
obvious (bearing in mind the response curve of the instrument) that longer periods 
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Fic. 7.—Variation of amplitudes and periods of short-period micro- 
seisms 1958 January 13 to March 7. 
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are associated with higher amplitudes. However, the true significance of this 
relationship is best shown by plotting individual amplitudes and periods directly 
against each other; this has been done in Figure 8. 

Three boundary phenomena emerge. 


1. There is a lower cut off period of 1s below which no activity of measurable 
amplitude occurs. This is merely the result of instrumental sensitivity. 
There is a higher cut off period of 3-5s above which no activity occurs in 
the short-period storms. 

There is an upper limit to the amplitude which is period-dependent, but 
activity of any period can occur at any level below this maximum. 
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Fic. 8.—Raw values of amplitudes and periods of short-period micro- 
seisms 1958 January—March. 


If the envelope of the points in Figure 8 is corrected for the characteristics of 
the seismograph the relationship between the period 7, and the maximum ampli- 
tude, Amax, which can occur at that period is obtained. This curve is shown in 
Figure g and states that Amax « 742. 

Previously Romney (1953) has examined the spectra of individual storms and 
has found that the microseismic spectrum at Pinewoods or Harvard is essentially 
continuous in the range from } to 5s, and that the amplitude increases with period 
at a rate proportional to the third or fourth power of the period depending to 
some extent on the meteorological conditions over the adjacent regions of the 
North Atlantic Ocean. For the storm that Romney quotes the power indices of 
the period are 3-9 at both Harvard and Pinewoods, values in accord with the 
fourth power relationship found in the present study. No attempt appears to 
have been made to extend Romney’s work. The concurrence of the fourth power 
relationship between amplitude and period in two different geographical areas 
suggests that the relationship may be of some importance, as any theoretical 
explanation of the generation and transmission of microseisms must account for 
this quantitative relationship between the two parameters. Alternatively it may 
be useful in regional studies of absorption of surface waves in rocks. 
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Fic. 9.—Variation of maximum ground displacement from 
rest position with period, short-period microseisms. 
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. 10.—Raw values of amplitudes and periods of long-period micro- 
seisms, 1958 March 20 to May 10. 
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Romney’s analysis relates to all periods occurring during a single storm where- 
as the relationship reported here concerns only periods and maximum amplitudes 
reported every six hours. However the maximum amplitudes on which the en- 
velope is based occurred during storms and furthermore the broad spectrum 
microseismic storms of the type shown in Figure 2(b) would show spectra very 
similar to those derived by Romney. 


The long-period microseisms. ‘The individual six-hourly readings of amplitude 
and period of the long-period microseisms have been treated in a manner similar 
to those of the short-period microseisms. In Figure 10 the raw data are plotted. 
The maximum amplitudes of the long-period microseisms increase with period. 
The envelope of the amplitudes has been corrected for instrumental response and 
the corrected variation of maximum amplitude, Amax, with period, 7, is shown 
in Figure 11. This curve states that Amax « T®®. Thus there is an interesting 
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Fic. 11.—Variation of maximum ground displacement 
from rest position with period, long-period microseisms. 


quantitative difference between the maximum amplitude variation of the short- 
period microseisms (Amax « T**) and the maximum amplitude variation of the 


long-period microseisms. The possible cause of this difference will be discussed 
in Section 8. 
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6. The weather situations which generate the microseisms 


The short-period microseisms. A comparison of the short-period microseismic 
activity with the surface weather charts prepared by the Weather Central shows 
quite clearly that the microseisms are produced by winds over the Ross Sea. These 
winds can be associated with local depressions in the Ross Sea or with the larger 
cyclones which circulate round the southern oceans just outside the periphery of 
the Antarctic continent. Figure 12 shows a typical surface weather chart of the 
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Fic. 12.—Surface weather chart of the southern oceans and Antarctica, 
12h G.M.T., 1958 January 29. 


southern oceans, drawn from observations at 12 hours G.M.T. on January 29, 
during the short-period microseismic storm of 1958 January 27-30. The de- 
pression centred at 63°S, 166° W, produced light south-easterly winds over the 
Ross Sea, but these are intensified by a subsidiary low over the centre of the Ross 
Sea. At Cape Hallett southerly winds varied between 20 and 4okt (10 to 
20m/s) for most of the duration of this storm. Again in mid-February a south- 
easterly air stream produced by a depression centred about 70°S, 155° W (Figure 
13) generated one of the major short-period storms of the summer (see Figure 7); 
southerly winds of about 30kt were recorded at Cape Hallett. A similar weather 
situation during 1958 March 4~5 produced the last major short-period storm 
before the winter freeze up. It will be noted that it is mainly south and south- 
easterly winds which are responsible for the major microseismic storms of the short- 
period type. The disposition and paths of the depressions are such that strong 
northerly winds are very rarely produced over the western part of the Ross Sea, 
or in the region of the continental slope. This fact is of some importance in con- 
nection with the study of the long-period microseisms. 

One feature of the short-period microseismic activity as yet unexplained is the 
rapid rise in activity of 5s period just before the major short-period storms. This 
phenomenon occurred before each of the three major storms of 1958 which are 
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Fic. 13.—Surface weather chart of the southern oceans and Antarctica, 
06h G.M.T., 1958 February 14. 
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Fic. 14.—Amplitudes and periods of the long-period microseisms, 1958 

January 24 to March 5. Q = quiet microseismic conditions, S = short- 
period activity masking long-period microseisms. 
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shown in Figure 7. Figure 14 reproduces the six-hourly amplitudes and periods 
of the long-period microseisms during January, February and early March 1958. 
The 5s activity intensified greatly during January 26, February 12 and March 3. 
On the day following each of these dates a major short-period microseismic storm 
began while the 5s activity was still increasing. Furthermore the §s activity does 
not increase in period with increase in amplitude which is unusual among the 
microseismic phenomena studied here. The usual behaviour is that shown by the 
parameters of the long-period storm of February 8-9 (Figure 14), and this is very 
similar to the behaviour of the short-period microseisms in Figure 7. 


The long-period microseisms.—These studies of microseisms are the result of a 
seismological observatory being established in a unique geographical situation 
during I.G.Y. Likewise as a result of I.G.Y. reliable maps of surface weather are 
available for Antarctica and the southern oceans which allow an attempt at the 
correlation of meteorological situations and the resultant microseisms. Cyclone 
paths around the Antarctic continent tend to follow mainly northwest to southeast 
trajectories. A study of the surface weather maps indicates that the long-period 
microseismic storms begin at Scott Base as depressions, below 970 mb at the centre, 
approach the N.W. area of the Ross Sea. The centres of the depressions are almost 
invariably between latitudes 62° and 65° as they pass longitude 165° E and precipi- 
tate the microseisms. Figure 15 shows the surface weather map at o6h G.M.T. 




















Fic. 15.—Surface weather chart of the southern oceans and Antarctica, 
06 h. GMT. 1958 March 14. 


on 1958 March 14 as the first big long-period storm of the season was approach- 
ing its maximum. A deep low pressure area (955 mb) is centred at 64°S, 175° E. 
Winds from this cyclone would create northerly swell moving into the Ross 
Sea at this stage. A subsidiary low pressure area in the Ross Sea immediately 
to the north of Scott Base produced short-period microseisms of about 3s period 
contemporaneously with the long-period storm. By ooh on March 16 the minor 
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395 
low had disappeared and the major depression had moved east of the Ross Sea 
and filled to g80mb; at this stage the microseisms had returned to background 
level. This situation was repeated during almost every long-period micro- 
seismic storm during 1958 March-May. Figure 16 shows the surface weather 
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Fic. 16.—Surface weather chart of the southern oceans and Antarctica, 


oohrs. G.M.T., 1958 April 20. 


situation at ooh on 1958 April 20 (storm of April 20-22) where the conditions 
are similar to those shown in Figure 15, the cyclone being centred at 61°-55 
166° E. 

The tracks of depressions of both the meridional and circular types follow a 
N.W.-S.E. course as they cross the Ross Sea (Alt, Astapenko & Ropar 1959). As 
a result of this, strong northerly winds are never observed in the N.W. part of the 


Ross Sea. Heavy northerly swell can only be generated to the north of the Ross 
Sea. 


7. The relationship of microseism and sea-wave periods 

The relation of sea-wave periods to periods of microseisms caused by these 
sca waves is of considerable importance. As is well known, the observation that the 
period of swell reaching the coast is double that of the prevailing microseisms 
(Deacon 1947, Bernard 1941) has been held to give strong support to the standing 
wave theory of microseisms. No wave measurements have been made during the 
1.G.Y. in the Ross Sea region but the parameters of swell will be discussed on the 
basis of the meteorological situation to see if the correlation of microseism period 
with either the period or the half-period of the significant sea waves leads to any 
absurdities. For this the theory of relations of wind, sea and swell developed by 
Sverdrup & Munk (1947) is used. 


The short-period microseisms.—The sharp upper limit of 35s to the period of 
microseisms caused by winds blowing over the Ross Sea must be due to some factor 
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associated with the sea waves such as limited fetch or maximum wind speed, or 
to a factor possibly connected with the second step of generation such as water 
depth. For deep water waves (where h > 4L, h being the depth of water and L 
the length of the wave) the period of the gravity waves is related to the velocity of 
progress of the wave by the expression t = 27C/g. A sea wave of period 3°5s thus 
moves at 5:5 m/s. According to Sverdrup & Munk significant waves of this period 


will be produced by 


(a) a wind of 10 m/s ( = 19kt) blowing for 6hr, or 
(b) a wind of 6m/s ( = 11 kt) blowing for 12hr. 


In the former case a fetch of 50km is required and in the latter a fetch of 
110 km. 

A sea wave of period 7-0s moves at 11-0 m/s and would be generated by a wind 
of 2om/s (38kt) blowing for 12br over a fetch of 180km. On a basis of un- 
limited duration a wind of 20 m/s would produce significant waves of 6-0, 7-4 and 
8-3s period at fetches of 100, 200 and 300km respectively. 

Thus in order to produce sea waves less than the maximum period of 3°5s in 
the Ross Sea the fetch must be less than 50 km, the wind speed less than 6-8 m/s or 
the duration of the storm less than 12hr. It is difficult to imagine that fetch, 
windspeed and duration of wind are limited to these maximum values during the 
major meteorological storms of the summer, which constitute the only occasions 
when the 3:5s microseisms were recorded. It can be seen (Figure 1) that an 
area of about 400km x 300km might come under the influence of a wind system. 
A reasonable fetch might be one-half of these dimensions, i.e. 150-z200km. The 
duration of build-up of the main microseismic storms is 24 hr. During the summer 
of 1957-58 mean daily windspeeds of 8—10 m/sec were usual at Scott Base during 
microseismic storms. As the low pressure centre was invariably some hundreds 
of miles to the north of the base on these occasions it must be expected that wind- 
speeds considerably higher occurred over the Ross Sea. 

Thus it seems virtually certain that in the Ross Sea waves of periods 
greater than 3-5s can be generated and that in fact 7s would be the best estimate 
from meteorological data of an upper limit during the summer storms. This means 
that the limit of 3-5s to the microseism period is due to a cut-off mechanism in 
the generation process, or to a subdivision of frequency in that process. 


The long-period microseisms.—The long-period microseisms have periods be- 
tween 4 and 1os. During summer and autumn the dominant periods are in the 
lower part of this band. ‘Table 1 gives the distribution of these periods between 
1958 January 1 and May 30. 


Table r 


Period 
Range (s) 


Frequency 


The long-period microseisms are believed to originate in the swell created by 
north to north-westerly winds of the cyclones approaching the Ross Sea from the 
N.W. As the centres of these cyclones pass 165° E microseisms of the long-period 
type appear at Scott Base. 
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According to Sverdrup & Munk waves of 8-gs period would be generated 
by winds of 20 m/s and waves of 10°6s by winds of 30 m/s, over a fetch of 
400km. The former velocity is equivalent to 8 Beaufort and the latter to 11 
Beaufort. 11 Beaufort is rarely recorded and then only for short periods so it is 
to be expected that waves produced by cyclones to the north of the Ross 
Sea will rarely exceed 10s in period. These periods are shorter than those re- 
quired to explain the microseism periods on the basis of the half-period relation- 
ship. However Darbyshire (1952) has also studied the generation of sea waves, 
by wind and, on a basis of empirical observations, gives periods to the waves of 
maximum amplitude produced by corresponding gradient winds which are higher 
than those of Sverdrup & Munk, particularly at the higher periods. According to 
Darbyshire’s data gradient winds between 30 and sokt (16 to 27 m/s) would pro- 
duce sea waves of which the period of maximum amplitude would vary between 
8 and 14s. This wave-period range encompasses the double period range of the 
summer long-period microseisms. 

Microseism periods were in the 7-10s range during the winter of 1957 not 
only at Scott Base but also at Hallett, Wilkes and Byrd Stations. The storms of 
the spring equinox 1957 can be correlated at all bases. ‘The maximum amplitudes 
were recorded at Wilkes Station which suggests that the origins of the microseisms 
were to the north or west of that station. The highest amplitudes of these winter 
storms at Wilkes Station (see Figure 19) were about 20~30 per cent of the ampli- 
tudes of the greatest microseismic storm at Wilkes when (on 1958 April 24-25) a 
depression was centred just to the north of Wilkes producing winds gusting to 
116kt at the Station. This 6~7s storm was hardly recorded at Scott Base, Hallett 
or Byrd Stations. This suggests that only the very long-period microseisms can 
effectively penetrate mountain ranges and other barriers so that such microseisms 
can be transmitted over much of Antarctica. The reason for this apparent region 
of generation being in the vicinity of Wilkes Station is discussed in Section 8. 


8. Discussion of the observations 


The primary source of microseisms.—During recent years most workers in the 
microseism field appear to have accepted the view that microseisms are caused by 
sea waves associated with cyclonic depressions, rather than true surf or atmospheric 
oscillations. There is, however, considerable opposition from a few workers to 
the uninhibited acceptance of sea waves as the universal primary source of micro- 
seisms. Donn & McGuinness (1960) correlate the passage of cold fronts seaward 
from the coast of North America with the abrupt commencement of microseisms 
at Palisades. They also show that high onshore progressive waves are not directly 
related to storm microseisms and in general negate any direct wave mechanism 
for the origin of the microseisms. They suggest that either direct coupling from 
the relatively high energy atmospheric turbulence or some form of resonant coupling 
from barometric oscillations travelling with acoustic speeds are more likely sources. 

The studies of microseismic activity at Scott Base are entirely in agreement with 
the concept of swell at sea as at least part of the primary source of microseisms. 
In particular it will be noted that, 


(a) Microseisms are only observed when cyclones, and thus swell producing 
conditions are present in the Ross Sea, or to the northwest of the Ross Sea. 
(b) If microseisms were due to atmospheric oscillations unassociated with sea 
roughness, they would not be materially reduced by an ice cover on the sea. But 
2C 








398 T. Hatherton 


the locally produced short-period activity is reduced to practically zero by freezing 
ef the Ross Sea in winter. 

(c) The maximum long-period activity occurs after freezing of the sea around 
the coast has begun, when conditions prohibit the production of true surf (and 
also the reflection of swell). 


Recently Santé (1959a) has examined a large number of Japanese records of 
microseisms and swells. He also concludes that “the direct cause of microseisms 
is undoubtedly swell.”” However, the evidence that Donn & McGuiness (1960) 
have provided requiring the coupling of atmospheric disturbances to the sea in 
the generation of microseisms is very strong. Consequently either two mechanisms 
exist for the generation of microseisms or alternatively both swell and atmospheric 
disturbances are required for the generation. Nanda (1959) has opened the dis- 
cussion of an “‘all embracing cause involving winds, waves and special orientation 
of winds...” 


The secondary source of the microseisms.—It is necessary to introduce the concept 
of a secondary source to describe the place where some of the energy of the sea 
waves is changed into microseisms. It has been suggested above that the long- 
period microseisms received at Scott Base are due to swell generated in cyclones 
to the northwest of the Ross Sea and which moves onto the continental area below 
the Ross Sea there producing the microseisms. This picture of events explains the 
discrepancy between the amplitude—period relationships of the short- and long- 
period microseisms. 

We have for the amplitude of the short-period microseisms, 


; o 7a" (1 ) 


and for the long-period microseisms 


Amax oc Tf. (2) 


If it is accepted that the long-period microseisms are created by swell approaching 
the continental shelf of the Ross Sea, the geometry of generation and decay being 
as shown in Figure 17, then the disparity between equations (1) and (2) can be 


62°s 


c's 7's 78s 
“ fetch - decay >be — Ross Sea —+ 
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Short- period 
microseisms 


Fic. 17.—Schematic illustration of the sea wave processes involved in 
the generation of the long- and short-period microseisms. 7, H indicate 
period and amplitude of long-period sea waves, t, h indicate period and 
amplitude of short-period sea waves, subscripts f and d indicate values of 
parameters at end of fetch and decay respectively. 
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explained in terms of frequency dependent decay of the sea waves responsible 
for the long-period microseisms. Little is known about the loss of energy of sea 
waves though most theories require that the decay be frequency dependent. 
Empirically Darbyshire (1957) finds that the decrease in energy is independent 
of wave period and that the broadening of the wave front is the most important 
factor contributing to this loss. Darbyshire’s data cover the wave periods 13 to 
18s which are mainly greater than the periods considered here. Nor is anything 
known about the role of the sub-zero air temperatures in the loss of energy of 
the waves. 
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Fic. 18.—Relative amplitudes of short- and long-period activity. 
(a) X = Amplitude of maximum 6s activity on short-period vertical 
component compared with extrapolated Amax-T curve of short-period 
activity. 

(b) X = Amplitude of maximum 3°5s activity on long-period N—S 
component compared with extrapolated Amax-T curve of long-period 
activity. 

@ Abcissal transposition of above points to compare amplitudes of hori- 

zontal and vertical motion. 


Figure 18 reproduces the maximum amplitude-period relationships of both 
types of microseisms. In addition it shows the maximum 6s activity recorded on 
the short-period vertical component at the height of the biggest long-period 
storm, and the maximum 3°5s activity on the N-S component at the height of the 
greatest short-period storm of 1958. Thus with the extrapolation of the Amax—T 
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curves the amplitudes of 3-5 and 6s activity produced by distant meteorological 
events can be compared with the activity produced by local events. Figure 18 
shows that the 6s activity produced by distant storms is about 30 per cent of 
the amplitude of the locally produced 6s activity if the extrapolation is valid. The 
distant 3°5s activity is only 14 per cent of the maximum local activity of that 
period. This would be satisfied by frequency dependent decay of the swell as it 
moves from the fetch to the Ross Sea. In tne case cited the heights of the waves 
at the end of the decay distance (Hq) relative to the heights at the end of the fetch 
(Hy) are proportional to about T;!%. This reconciles, approximately, equations 
(1) and (2). 

Abscissal transposition in Figure 18 of the amplitudes of the 3-5 and 6s activity 
from the component on which these amplitudes were measured to the other com- 
ponent indicates that the horizontal component of motion in the N-S plane is 
about double the vertical component. 


Transmission phenomena, microseismic barriers and the continental shelf.—It 
appears that the data so far described are entirely in accord with the idea that 
microseisms of large amplitudes can only be generated on the continental shelf. 
There is a growing body of opinion that this concept is in general correct. The 
two most recent major summaries of large-scale microseismic investigations, in 
America and in Japan, state that, “for a continental station it (the generating area) 
is usually within 200 miles of the shore” (Carder & Eppley 1959, p. 186), and “the 
energy transformation from swells to microseisms occurs somewhere near the 
coast” (Santo 1959b, p. 494). These conclusions were anticipated by van Straten 
(1953) on the basis of a review of evidence available at that time. 

The hypothesis of the role of the continental shelf was tested by carrying 
out seasonal studies of microseismic activity of the long-period type at Wilkes 
(66° 15'S, 110° 30’ E) and Hallett Stations respectively. These bases lie just on the 
continental side of the shelf and the seasonal control of the sea surface by the pack 
ice should give seasonal variations in microseismic activity similar to that at Scott 
Base. Figure 19 shows the daily mean amplitudes of the long-period microseismic 
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‘IG. 19.—Seasonal variation of long-period microseismic amplitudes, 
(a) Hallett Station, (b) Wilkes Station. 


activity at Wilkes Station and Hallett Station from 1957 August 1-1958 July 31, 
and these can be compared with the long-period Scott Base results for the same 
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year. Hallett Station shows exactly the same seasonal variation as Scott Base and 
the microseismic storms are identical in times of occurrence and in relative ampli- 
tudes. Figure 20 compares the periods at Scott and Hallett at the same time and 
indicates that the periods are similar and justifies the period selection procedure, 
for the periods were read not only from widely different instruments but also 
from different components. The Scott periods are about 0-2s less than the Hallett 
periods at 8s but about o-2s longer at 5-5s. 
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Fic. 20.—Differences between mean daily microseismic periods at Scott 
Base and Hallett Station. 


On the other hand, the Wilkes Station variation shows marked differences 
from those of Scott and Hallett. Taking the amplitude of the maximum storm at 
each of the three stations as unity to overcome problems of instrumental amplifica- 
tion it can be seen that storms of large amplitudes occurred at Wilkes during August 
and September 1957, when the edge of the sea ice theoretically was at its northern- 
most limit. Periods during these storms were of 73-108; the microseisms were 
also recorded, though of lower amplitude at Scott and Hallett. The microseismic 
storms of March to May are of shorter predominant period and a frequency analysis 
is required to determine whether the spring microseisms also contain a substantial 
74-10s component. These autumn microseismic storms at Wilkes cannot be cor- 
related with the autumn storms at Scott or Hallett. The equinoctial appearance 
of the spring microseisms is probably related to the observation due to Schwerdt- 
feger (M. Rubin, personal communication) that the meridional temperature 
gradients and, hence, the zonal westerlies appear to have a maximum at the time 
of the equinoxes. 

Surface waves decreasing in period from 20 to 7}$s are produced at Wilkes 
and Hallett Stations and at Scott Base by earthquakes with epicentres in the eastern 
part of the Indian Antarctic swell. A long coda of 74s waves with strong vertical 
motion is also a characteristic feature of these shocks. These waves have not yet 
been analysed but Ewing & Press (1956) have ascribed the short-period charac- 
teristics of other examples of this infrequent type of seismogram to short-period 
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Love waves; scattering and transformation of the Love waves in the vicinity of 
the seismograph station is held to account for the vertical and longitudinal com- 
ponents of motion observed on this type of seismogram. Whatever the detailed 
nature of these waves the presence of periods down to 7}s makes it feasible that 
microseisms of periods greater than 74s may be produced at storm centres in the 
S.E. Indian Ocean and transmitted via oceanic crust to the Antarctic Continent. 
Strong microseisms of periods less than 74s may only be produced by ocean swell 
over the continental crust. 

Iyer and others (1958) strongly deprecate the requirement of a “‘microseismic 
barrier”, preferring to explain anomalies and regional differences in the reception 
of microseisms in terms of refraction effects produced by sea bottom topography 
on Rayleigh waves (called the Stoneley effect by these authors). However, in the 
case of long-period microseisms at Scott and Hallett there appears to be a very 
real increase in microseismic activity when removal of the ice cover allows access 
of storm or swell to the continental shelf or slope. 


The amplitude-period relationship and the standing wave theory.—It has been 
shown that 


Amax « T4 (3) 


for the short-period microseisms; this relationship is also true for the long-period 
microseisms if these microseisms are due to swell moving onto the continental 
shelf from the north. 

The standing wave theory of microseisms developed by Longuet-Higgins 
(1950) is the most satisfactory, and indeed the only quantitative theory yet derived 
which gives a convincing physical mechanism for microseismic generation. This 
theory which explains the half-period relationship between the sea waves and 
microseisms, and which accounts for the amplitudes of ground motion observed 
is widely accepted particularly by workers in Great Britain. However there is still 
opposition to the theory by many seismologists, this opposition being mainly due 
to a reluctance to accept standing waves at sea as a common occurrence. For 
instance, it is very difficult to envisage any process which derives standing waves 
from progressive northerly swell entering the Ross Sea in April, when the coast- 
line is protected by sea ice. Consequently further tests of this theory are desirable 
if it is possible to devise them. Longuet—Higgins has shown that the mean pressure 
p: on a plane horizontal sea floor at depth z below the undisturbed ocean surface 
for two oppositely directed waves of equal frequency is given by 


Pz = Pst gpz—dhyheo® cos 20t 


where ps denotes the pressure at the free surface (assumed constant) p is the 
density of the ocean, /), ho are the amplitudes of two sinusoidal waves of period 
2n/o. 

Longuet- Higgins has also shown that the more widely the energy is distributed 
in the spectrum the smaller is the resulting disturbance. ‘The maximum pressure 
variation on the sea-floor would occur when two oppositely directed wave trains 
of single frequency meet with amplitudes equal and at a maximum for the period 
considered. Neumann (1952) has given data relating maximum amplitudes and 
periods of fully developed sea waves for an unlimited fetch, and shows that h oc ty1°7. 
More recently Darbyshire (1959) has found that h, ot,3 where t, and h, (some- 
times called the significant wave period and height) refer to the mean period and 
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height of the highest third of the waves. Hence using the latter relationship, 
Pz—Ps 


ee (4) 


or the maximum amplitude of the microseisms on generation «T?. 
The translation of (4) into (3) depends on several factors, principally 


(i) the period dependence of the size of the generating area, 
(ii) the propagation characteristics of the microseisms, and 
(iii) the attenuation of the microseismic energy. 


None of these factors can be discussed quantitatively with any confidence. If 
it is assumed that the generating areas are inversely proportional to wave number 
and that the propagation characteristics are independent of frequency, attenuation 
must be responsible for the translation of (4) into (3). If the attenuation coefficient 
is of the form k = 7/OTV, where V is the velocity of the microseisms (= 3-6 km/s), 
the value of O which most satisfactorily explains the frequency dependence of the 
amplitude decay is ~ 50. This is very low compared with the commonly quoted 
value for very long waves of 150 (Ewing & Press 1954). Darbyshire (1957) states 
that the period of a sea wave depends on the square root of the wind speed and that 
the amplitude of the wave depends on the square of the wind speed. A combination 
of these relationships gives h o ty’, which makes the fourth-power microseismic 
equation compatible with the standing wave hypothesis and a O of 150. Un- 
fortunately the unsatisfactory state of quantitative knowledge of all the above 
factors hardly makes a critical test possible at present. 

Nevertheless the fourth-power relationship between microseismic amplitude 
and period may be of use when other factors can be properly assessed. Meanwhile 
it will be interesting to see if such a relationship is general, particularly at shore 
stations where the high-frequency microseisms are still observable. 

The role of the water depth.—Several authors have discussed the effect of the 
depth of water (a) under the generating area and (b) between the generating area 
and the station. It follows that if the microseisms are Stoneley waves (Stoneley 
1926) the water depth has a considerable effect upon their transmission. Mitra 
(1957) has investigated the propagation of Rayleigh waves in a multilayered medium 
with the uppermost layer of water. Associating the prevailing periods of micro- 
seisms with stationary values of the group velocity Mitra finds that the predominant 
periods depend critically on the nature of the first layer of the basement, but at 
depths large compared to the thickness of the first layer both group velocity and 
period are practically independent of the nature of the layer. Imbert (1954) has 
demonstrated a linear relationship between the period and the depth of water 
(plus unconsolidated sediments) in the generating area. 

If the physical picture of microseism generation in the Ross Sea is that hypo- 
thesized in this paper it is difficult to ascribe any significance to the depth of water 
except in its general relationship to the continental shelf. The maximum amplitude 
has been shown to vary in a continuous manner with period which is consistent 
with the increase in height of ocean waves with period, accentuated by frequency- 
selective generation and absorption processes. Further single storms almost 
invariably show an increase in period with increase in amplitude as a cyclone 
deepens. There are, however, several factors which may be associated with water- 
depth control. Among these are (a) the 5s microseismic storms that precede the 
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major short-period activity (referred to in Section 6) and (b) the apparent sharp 
cut-off activity above 3-5s for the short-period microseisms and above 7s for the 
autumn long-period microseisms. The latter does not seem in accord with the 
observed decline of low frequencies in sea wave spectra. 

To investigate these points and to clarify other issues discussed in this paper 
a more sophisticated frequency analysis is necessary. 


g. Conclusions 


1. The seasonal variation of ice cover over the Ross Sea and the oceans sur- 
rounding the Antarctic continent offers a partial control over the complex wind 
sea—crust phenomena which generate and transmit microseisms, and this may be 
useful in the investigation of both microseisms and crustal structure. 

2. Three classes of microseisms are recorded at Scott Base: short-period 
microseisms (1-3°5s) which are due to winds over the Ross Sea; long-period 
microseisms (4~7s) which, it is thought, are created by swell from meridional and 
circular cyclones moving onto the continental shelf in the Ross Sea; and ultra 
long-period microseisms (7$—10s) which are considered to be produced at or near 
the storm centre and transmitted through both oceanic and continental crusts to 
the Base. 

3. The maximum amplitude of ground motion varies as the fourth power of 
the period of the microseisms. This relationship may allow a test of any quantitative 
theory of microseisms though at present the uncertainty in the mode of transmission 
and the lack of reliable attenuation measurements make this test difficult to apply. 
Confirmation of a mechanism of generation would, conversely, allow an investiga- 
tion of transmission phenomena. 

4. The complete absence of short-period microseisms when the Ross Sea is 
frozen confirms that sea roughness is essential to the generation of microseisms. 
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Reflexion and Refraction of Plane Elastic Waves at a 
Plane Boundary between Aeolotropic Media 


M. J. P. Musgrave 


(Received 1960 July 19) 


Summary 


General equations in stress and displacement are set down and 
their implications qualitatively discussed. 

Detailed results are presented for the slownesses, amplitudes and 
energy fluxes of body waves generated by the incidence of a body wave 
upon a range of differently oriented boundaries in hexagonal media. 

Results for ice and beryl show that assumption of isotropy preserves 
the qualitative form of the reflexion characteristics; in contrast, the 
variations caused by changing the orientation of the boundary in zinc 
are too great to be adequately represented in terms of isotropic constants. 


1. Introduction 

In recent papers (Musgrave 1954 a, b, Synge 1956, 1957 and Buchwald 1959) 
the dynamics of elastic waves in aeolotropic media has been presented from 
various aspects and quantitative examples in terms of hexagonal media have been 
given. The general reflexion-refraction problem for plane elastic waves at a plane 
interface between aeolotropic media has been sketched in recent reviews (Musgrave 
1959, 1960) but no quantitative information appears yet to have been made avail- 
able. The purpose of this communication is to present a qualitative discussion of 
the general problem together with some numerical results obtained for the reflexion 
of body waves at the free surface of some hexagonal media. 


2. Notation 


For the general problem we shall use notation derived from Synge (1956); for 
the particular case of hexagonal media, a non-tensorial notation is more conve- 
nient. Following Synge (1956), we may represent a plane progressive wave by the 
real or imaginary part of the displacement 

Uz = UR*+iugt = Ap, exp tw(sjxj— 2), (1) 
where px is the unit displacement vector associated with the slowness sz. Substitu- 
tion of the form (1) in the equations of motion 

9ij,f = CijklUn,1j = Pui, (2) 
yields 


[pdix— 4x ]pi = 0, (5%, Kronecker delta) 
406 
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where aj = $45:Cizx1. Thus sy and px are obtained as characteristic roots and 
vectors which may be real or complex. If px and sy are real, ux* or ux’ represents 
a wave of linear polarization and zero attenuation; if px and sz are complex, u,* 
or u,' represents an elliptically polarized wave whose amplitude is attenuated with 
distance. 
The progressive flux of energy F;* which is a time average of — o,;*u;* is given 
by 
Fy* = $A? wei sual pjsxpit Pisepr), (4) 


where the bar indicates a complex conjugate. 


3- General conditions for displacement and stress at a plane 
interface 


Suppose two semi-infinite aeolotropic media are in contact in the plane xg = o. 
Let their densities be p and p’ and their Hookean elastic constants referred to the 
rectangular co-ordinate systera x; be cjg: for x3 > 0 and c’gjx; for x3 < 0 
respectively. 

If we represent an incident wave of unit amplitude in the first medium by a 
displacement of the form 


uc) = pe exp iw(sjx;— 2), 

the possible reflected and refracted waves will have displacements 
Ap, exp iw(sjx;—t), 

and 


A’) py’ exp iw(s;'x;— t), 
respectively. 

At this stage, we have not disposed the x; and x2 axes uniquely. In formal alge- 
braic terms, since the c;jx; represent twenty-one constants, zero or non-zero in 
magnitude, no generality is lost by choosing the x,-axis perpendicular to the 
normal to the boundary x3 and to the incident slowness vector; the x-axis is then 
the intersection of the plane of incidence and the boundary. It should be pointed 
out that this choice of axis is not necessarily the most convenient if quantitative 
calculations are to be made. 

For a boundary at which no slip occurs, we require continuity of displacement 
and stress over x3 = 0. In order that these conditions are obtained for all times, 
the phase of all the waves at the boundary must be equal at any instant ft; thus 


so = sg) = $9') (5) 


and these relations in fact determine the number of possible reflected and refracted 
waves. Consider Figure 1: TO is the trace on the boundary of the incident slow- 
ness. When the energy associated with the incident mode meets the interface, a 
redistribution of the modes of propagation takes place in accordance with con- 
tinuity of displacement and of stress. Conditions (5) require that the traces of the 
slownesses of possible reflected and refracted waves are all equal so that the 
periodicity at the interface is undisturbed. Hence OT’ represents the common 
trace of all the possible reflected and refracted waves and their slowness vectors 
have extremities defined by the intersections, real or imaginary, of the normal 
to the interface through 7’ with the slowness surfaces of each medium. These 
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slowness vectors have the form (0, se, s3(4”)) and (0, so), s3'(%)) where M, N 
= I, 2,...6since the slowness surfaces are of the sixth degree. [The familiar fact 
that the planes of incidence, reflexion and refraction are coincident (i.e. x1; = 0) 
is implicit in these forms. However it is essential to appreciate that these planes 
are defined by the wave normals or slowness vectors and not the rays or energy 
fluxes about which entities no such simple generalization is possible in the case of 
aeolotropic media. Incident reflected and refracted rays are, in general, not co- 
planar.] For s3™), s3’(%) real, the reflected and refracted waves are body waves 
and the centro-symmetric, continuous character of the three sheets of the slowness 
surfaces results in the normals at three points of intersection on each surface having 
a positive x3-component, say M = 1, 2, 3, N = 4, 5, 6, and the remainder have a 
negative xg-component, say M = 4, 5, 6, N = 1, 2, 3. Hence there are three 
possible waves with appropriate slowness vectors and energy fluxes for each 
medium. 





Fic. 1.—The construction in slowness space yielding the extremities 
of the slowness vectors of the three reflected and three refracted waves. 


When s3(™), s3’(") are complex, they must occur in conjugate pairs say M, N 
= 1 and 4, 2 and 5, 3 and 6. These waves have normals of type (0 : sg) : s3*(4) 
and their amplitudes suffer exponential attenuation or amplification exp( —wsg'x3(™)) 
as well as sinusoidal variation exp(iwsgx3) with distance from the boundary (Stoneley 
1955, Synge 1957). Accordingly, from any such pair the wave representing 
attenuation of the amplitude is retained. The energy flux for such waves has been 
shown to be parallel to the boundary by Synge (1956). 

There are thus three slownesses in each medium which provide for physically 
realizable waves and whose amplitudes remain to be determined by the boundary 
conditions in displacement 


3 3 
PeP+ > AM) pM) _ > A’ py"), (6) 


M=1 N=1 
and in stress 


3 3 
( WN 
o3f2+ > 03 = > oj, 
M=1 N=1 
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where 
03 ju >= twAMc, peipe™! 5) exp tw(S¢ M ‘Xq—t ) 
and ten 
037) = iwA'Ne'g ip’ 9’ exp iw(sq'xq— 1). 


The six equations, in general complex, represented by (6) and (7) uniquely 
determine the amplitudes A(”), A’(%) (M, N = 1, 2, 3). Solutions for these 
amplitudes may be checked by using the balance in energy flux normal to the 
boundary 

3 3 
F3¥(D + * F3*) ; F3/*), (8) 


—_— 


M=1 N-1 


A balance in energy flux parallel to the boundary may also be considered. However, 
when due allowance is made for different cross sections through which beams of 
energy having a common area at the boundary must pass, equation (8) is again 
obtained. 

It will be apparent that the general problem is fully soluble though somewhat 
tedious. In general, there is a choice of three types of incident wave whose angle 
of incidence may be varied and for each such set of incidence conditions the 
orientation of the medium on either side of the interface may be varied. Clearly, 
for any given material the variation of these parameters alone could yield a set of 
data so vast as to be quite inappreciable and merely embarrassing. 

For that reason, in order to obtain a small insight into the effects of anisotropy 
on waves incident upon a boundary, a much simplified problem has been examined: 
the reflexion of waves at the plane surface of a hexagonal medium when the 
plane of incidence contains the zonal axis. 


4- Reflexion of waves in hexagonal media 


We consider a plane boundary of arbitrary orientation in a hexagonal medium 
and let the plane containing the normal to the interface and the zonal axis be the 
plane of incidence. 

Let the space co-ordinates x; retain their significance and let the coordinate 
system (x, y, z) be referred to the zonal axis and two directions in the isotropic 
plane. For the problem in hand we may let x = x; and the direction cosines of 
0x3 with respect to oy, oz be (0, m, n), where n = cos @ as shown in Figure 2. 
By analogy with previous notation referring to the axis of symmetry in hexagonal 
crystals, we denote the displacement vector (p(™%), g(%, r(%) and slowness 
(sz, sy, s(™)) where M may designate the quasi-longitudinal (L) or quasi- 
transverse (7), T2) character of the wave. 

We may now express the stress conditions at the boundary in terms of quantities 
already calculated in terms of the five independent elastic constants of a hexagonal 
crystal. 

Thus the slowness components are related by 


S2 = NSy— MSz, §3 = MS y + NSz 


and the non-periodic magnitudes of the stresses are given by 


o13™) = mx + nxz™, 


cog”) = m n(yy) —_ 22M) — (m2 — n?)y2) 


a3g™) - m2yyM) + n2z2) + 2mn yr), 
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where 


yy) — BOO(61 qs, + c1g7)s,(M), 
zo) — BW(ey5qs, + cggr'Ms 0), 

VE) = BWMeq4(1Ms (MD 4 GMs), 

3M) = BWM eq4( PMs (M) + Ms (M)), 

xy = BW3(ey) a €12)(q™s2 + ps), | 

BMD — iwA™), 
the appropriate quantities being used to calculate the stresses due to the wave 
designated by (M). Now consider an incident wave of type 7) which in hexagonal 
media possesses a truly transverse displacement (1, 0, 0) in the present co-ordinate 


systems. This wave generates only the stress o3 at x3 = o and hence only a 
single reflected wave of the same type 7; where 


so) = 5(R) 


and s3‘®) is the appropriate root of S(7)(soW), s3(*)) = o, which represents the 
section in the plane of incidence of the spheroidal sheet of the slowness surface. 


Xo, m,n) 


Incident 





Fic. 2.—The disposition of reference axes in the problem of reflexion at 
the boundary of a hexagonal crystal. 


The equation 
ais) + aig = 0, 


then yields the amplitude of the reflected wave as 


A(R) — m(C11 — C12)Sy + 2ncagsz? 


m(C11 — C12)Sy + 2ncags es” 


Since (sy), sz“) and (sy‘*), s2(®)) both lie on a zonal section of the separable ellip- 
soidal sheet, they each satisfy the equation 


(€11 — €12)Sy? + 2C448,7 = 2p. (14) 
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In conjunction with the condition s2¥) = s9(*), which may be written 


n(sy) —sy(®) = m(s) —s,(P), 


it may be shown that 
m(C11 — €12)(Sy + sy”) + 2ncgy(sX +5”) = 0, (15) 
and hence A® = —1 for all,incident waves. 

Two other types of incident wave are possible, L and Tz. Each of these gives 
rise to stresses o23 and o33 at x3 = o and in general results in reflected waves whose 
slownesses are determined by the phase condition at the boundary and whose 
amplitudes are given by the equations in stress 


025 + ong + ao37) = 0, | 


(16) 


o33 + 033) + o337» = ©, | 


Amplitudes so calculated must also satisfy the equation in energy flux 


Fx? + FX + FyT) = 0. (17) 


Numerical calculations of the reflected body waves due to a wave (L or 72) of 
unit amplitude incident upon an arbitrarily oriented plane boundary in zinc, ice 
and beryl were made. With the exception of the graphical solution of the phase 
condition, the whole calculation was programmed for the peuce. A preliminary 
programme provided the data for plotting the slowness sheets. Angles of reflexion 
for given incident waves were then found graphically. A main programme then 
provided for the calculation of stresses, the solution of the equations for the 
amplitudes, the calculation of the components of energy flux Fz and F3. These 
last were used as a check on the basis of equation (8). Slight discrepancies in the 
check have been shown to be due to the relatively low accuracy of the graphical 
stage in the calculation. 

The use of a graphical method of solution for the possible reflected waves pre- 
cludes the possibility of finding complex roots for the slowness vectors which 
correspond with the surface waves generated by waves incident at angles greater 
than the critical angle. ‘To solve for complex slowness vectors would require a 
more elaborate DEUCE programme and this at present has not been urged though it 
is possible that the phase changes which occur when the critical — cle is exceeded 
may be of siguificant interest. 

A selection of results relating to the reflexion of body waves at a free boundary 
is here presented graphically. Figure 3 shows the amplitudes of reflected waves 
at a plane boundary in an ice crystal for (a) incident L-waves, (b) incident T-waves 
of unit amplitude; Figures 4(a) and (b), 5(a) and (b) show similar curves for bery] 
and zinc. The dotted curves in Figures 3 and 4 and the curves of Figure 5(c) show 
equivalent results for some hypothetic elastically isotropic media whose elastic 
constants are equal to certain of the crystal constants. Thus in Figures 3, 4 and 
5(c) curves I and IL were obtained by setting cy) = c33 = ¢1)'%, cag = gq, 
C12 = C13 = C11 — 2€44 = C4 49 — 204449 for I and similarly but with C11 = C323 
= ¢33'*9 for II. Curve III in Figure 5(c) is obtained by setting cy2 = ¢y3 = cy2(*, 
C44 = c4g(t and C11 = C33 = C12 + 2C44 = Cy9'X) + 2c44(X and IV similarly but with 
C12 = €13 = €13'%9, Figure 6 shows the relative magnitudes the energy flux normal 
to the boundary for reflected waves in zinc. The incident energy flux is here nor- 
malized to unity and the sum of ordinates for any given angle of incidence should 
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Fic. 3.—Amplitudes of waves reflected from boundaries of various 

orientations for (a) incident L-wave, (b) incident T2-wave, in ice; the 

dotted curves show amplitudes obtained in hypothetical isotropic media 
having some properties in common with ice. 
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Fic. 4.—Amplitudes of waves reflected from boundaries of various 

orientations for (a) incident L-wave, (6) incident T2-wave, in beryl; the 

dotted curves show amplitudes obtained in hypothetical isotropic media 
having some properties in common with beryl. 
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Fic. 5.—Amplitudes of waves reflected from boundaries of various 
orientations for (a) incident L-wave, (6) incident T2-wave, in zinc. 
(Fig. 5c) shows amplitudes obtained in hypothetical isotropic media 


having some properties in common with zinc. 
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Fic. 6.—Relative magnitudes of energy flux normal to boundary for 
reflected waves in a zinc crystal. 
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Fic. 7.—Unusual disposition of incident and reflected rays which can 
occur upon reflexion at a boundary in a zinc crystal. 


therefore be unity. Figure 7 shows a somewhat unusual disposition of incident and 
reflected rays which can occur at a boundary in a zinc crystal. ‘Table I gives the 
physical constants of the crystals used to provide the results from which the graphs 
have been prepared. The relevant sections of the slowness surfaces of ice, beryl 
and zinc are shown for reference in Figure 8. 


Table 1x 


Crystal cll C33 c12 C13 C44 


Ice 1°385 1°499 0°707 0°581 0°319 
Beryl 26°94 23°63 9°61 6°61 6°53 
Zinc 14°30 5‘0 1°70 3°30 4°0 


All ci are in units 10! dyne/cm?; p in g/cc. 
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Z-axis 








Fic. 8.—Sections in the zonal plane of the slowness surfaces of (a) ice, 
(b) beryl, (c) zinc. Units are 10-6 cm—! sec. 





418 M. J. P. Musgrave 
5- Discussion 


This problem was undertaken in order to investigate the effects of anisotropy 
in a simple reflexion problem. Deviations from results obtained on the assumption 
of isotropy are of interest and could be of significance to seismologists and users 
of ultrasonic methods. 

In the case of only slightly anisotropic media, such as ice and bery], it is clear 
that the assumption of an equivalent isotropic medium may offer a reasonable 
approximation to the actual reflexion characteristics of these materials although 
significant quantitative errors arise for certain ranges of incidence, and in relation 
to precise seismological measurements in regions where elastic anisotropy is appre- 
ciable, the stricter theory may be required. On the other hand, it would appear 
that the reflexion characteristics of a crystal like zinc vary so considerably with the 
orientation of the boundary, that no simple approximation can be regarded as 
acceptable even for qualitative purposes. 
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Summary 


In all seismic refraction surveys, the problem is to determine the 
constants in a system of equations of the type 


tii = aji+b;+ Ai; v 


where a; and 6; are ‘‘time terms” which are characteristic of the shot- 
point and seismograph station respectively, Aj; is the distance between 
the shot point and the seismograph, t;; is the time of propagation of a 
refracted wave and v is the velocity of propagation of seismic waves in 
an underlying marker layer. It is shown that the equations can be 
solved for interpenetrating networks of shot-points and seismographs 
provided that certain general conditions are satisfied. Factors which 
determine the uncertainties of the final solution are discussed, and 


methods of correcting for the effects of steeply dipping boundaries are 
included. 


1. Introduction 


During the last few years, the refraction surveys of the Dominion Observatory 
have been planned to an increasing extent with the idea of determining the constant 
‘time terms’ which enter into the travel-time equations for refracted waves. 
Examples of these projects have been described in various papers (Willmore & 
Scheidegger 1956; Millman, Liberty, Clark, Willmore & Innes 1960, etc.). A 
theoretical paper, describing a method of solving the equations has also been pub- 
lished (Scheidegger & Willmore 1957). This will be referred to later as Paper I. 

The essential feature of the method is that the time t;; which is required for a 


refracted wave to pass between two points of the Earth’s surface may be written 
in the form 


tig = aat+bj+Agy/v (1) 
where a and 5; are “time terms” which are characteristic of the shot-point and 
the seismograph station respectively, Ai; is the distance between the shot- point 
and the seismograph, tj; is the time of propagation of a refracted wave, and v is 


the velocity of propagation in a uniform underlying ‘‘marker layer”. If certain 
fairly broad conditions are satisfied, it is possible to determine the values of a;, 


* Published by permission of the Deputy Minister, Department of Mines and Technical Sur- 
veys, Ottawa, Canada. 
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b; and v, which gives the best fit to the observed travel-times tj. The advantages 
of the method are that the equations can be solved without requiring the shot- 
points and stations to be laid out in any particular pattern, that the maximum 
amount of information is extracted from the data, and that the requirements for 
making simplifying assumptions about the geological structure (assumptions which 
are, for example, involved in describing the structure in terms of plane layers) are 
minimized. 

In the course of various informal discussions, it has become evident that the 
method requires further clarification. In particular, we should like to demonstrate 
that it does represent a significant departure from conventional refraction methods, 
even though it uses the same basic equation. We shall investigate the extent to 
which errors in the observations or departures from the underlying assumptions 
are revealed by inconsistencies in the solutions, and we shall draw up a set of 
working rules which should be followed in all types of refraction experiments if 
uncertainties in the final solutions are to be minimized. With these ends in view, 
we shall first present the theory of the method in a more elementary form than 
that given in Paper I, and will then compare our recommended procedures with 
others which have been discussed in the literature. 


2. Solution of the time-term equations 


Following the procedure of Paper I, we shall assume that m recording stations 
observe waves from a total of m sources, so that up to mm equations of type (1) will 
be obtained. 

Apart from v, the problem contains (n+m) unknowns. These unknowns can- 
not all be determined uniquely, for if an arbitrary constant « is subtracted from all 
the’shot terms a;, and added to all the station terms 5;, the observational equations 
are unaltered. 

The observational equations are reduced to a set of (m+m) normal equations, 
each of which is obtained by adding together all the observational equations which 
contain observations from a given station, or of a given shot. The normal equation 
referring to the ith shot may be written 


jel jel 
mi(ai—%) = > (tis—Aj/v)— > (63+), (2) 


j=m j-m 


where m; is the number of stations which successfully observe the ith shot, and 
the summation signs on the right-hand side include the data from these successful 
stations. Similarly, the normal equations for the jth station may be written. 


i=1 i=1 


nj(bj+x) = > (tis—Ay/v)— > (i —2), (3) 
i=n i=n 


where nj; is the number of shots which are recorded by the station. Equations (2) 
and (3) together are equivalent to equation (2.3) of Paper I. 

In general, the solution of the normal equations requires the inversion of a 
matrix or of a cracovian, which may be prohibitively laborious if an electronic 
calculator is not available. The procedure, however, is greatly simplified if all 
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stations observe all shots, for in this case m = m and nj = n. We can now solve 
for a; by eliminating b; between (2) and (3) as follows: 


f-1 j=l i=l 


a me ES usu) >, (tu Sule i 52, a0) (4) 


=m j=-m imn 


In view of the indeterminancy of the system of equations, we may assign one 
time term arbitrarily, and we do this in such a way as to reduce £(aj—«) to zero. 
Equation (4) then becomes 


j=l i= 1 


j=l 
a= =e (tij — Agy/e )-—> m3 (tig — Aiy/v) + (5) 
j=m 


j=m i-n 
and equation (3) becomes 


{=1 


I 
b= - tes — Ays/v) —a&. 
ee” > (ti ij/V) —o 


i=-n 


On entering the time terms into the observational equations, we obtain nm 
equations of the form 


=1 j=l t-l =1 


LS (tiy-Ayle)-— > > (te — Agy/v) + =>, u-duio 


j=m j=m i=n 
= ty—Ay/vo—8y, (7) 
where 64; is the residual. This is an equation for 1/v which may be put in the form 
ciy—(1/v) diy = 845 (8) 


which is the same as equation (3.13) of Paper I. As a running check on the calcula- 
tion, we note that the sum of the residuals must vanish for each shot and for each 
station for all values of v, and hence that the sums of cj; and dj; must vanish 
separately. It also becomes evident that the cracovian or matrix, which is needed 
in the general solution, is simply a device for assigning the proper weights to the 
connections between the survey points when some of the possible observations 
are missing. 

If a few observations are missing from a substantial survey, it may be sufficiently 
accurate to determine v and most of the time terms by applying the short method 
to the survey points for which complete observations are available, and then to 
determine the remaining time terms directly from the observational equations. 

By the method of least squares, the best solution for v is 


I > Cij di; ( ) 
Uv io 3 a2; . 9 
On entering the value of v in equations (5) and (6), we determine all the time 


terms subject to the ambiguity which arises from the existence of the arbitrary 
term «. 
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In determining «, the best procedure is to arrange for some of the survey 
points to be used as both seismograph locations and shot points. This will yield 
a number of cases for which a; and 5; refer to the same point on the ground, for 
each of which we determine «jj; by setting a; = 5; in equations (5) and (6). The 
mean of all the available values of aj; is taken to be the best estimate of « for the 
entire survey. In surveys for which this procedure cannot be followed, the arbi- 
trary constant may be chosen to fit the known geology at some point, or to minimize 
implausible breaks in the time-term pattern. 

In principle, the time-term method can be extended to allow for a curvature 
of the Earth, or for an increase of velocity with increasing depth in the marker 
layer. In this case, terms of the type A;j/v will be replaced by a polynomial in 
Ajj, and the terms of this polynomial will be summed in forming dj;;._ It will still 
be possible to determine the coefficients of the polynomial from equation (8), al- 
though the solution will, of course, be more cumbersome than in the linear case. 


3- Corrections for inclined ray paths 


The foregoing theory has assumed that a satisfactory approximation for the 
travel time ¢ is (from Paper I): 


A H,acos@ Hpycosé 
t=—+ + 
v2 U1, V1 





’ (10) 


where @ = sin~1(v/v2), v1, v2 are the velocities above and below the marker, and 
Ha, Hp are the lengths of the perpendiculars drawn to the marker from the shot 
point and seismograph respectively. 




















Fic. 1.—Ray paths through a structure with steeply dipping sides. 


Let us now consider the more general case shown in Figure 1, in which a 
seismic ray originates at A, enters the marker layer at C, emerges from the marker 
at D and is recorded at B. CE and DF are tangents to the marker, AE and BF are 
normals to CE and DF, and G and H are the feet of the normals from C and D 
to EF. Let AB = A, EF = A’, AE = Hy’, BF = Hp’. Hg’ and Hy’ will be close 
approximations to H, and Hp respectively, unless the interface is very strongly 
curved in the vicinity of C and D. 

It will now be shown that a much closer approximation to the travel time is 


Hi cos Pia Hy’ cos 6; 
~ . 





" (11) 


Vi 


where the angles @iq and 4)» are shown in Figure 1. 
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The proof involves only the assumption that CD is parallel to EF. Then 


AC CD DB 
itera Amini 

V1 v2 vi 

A’ AC EG DB HF 
+(— om a 


‘9 


V2 V1 vd 

, H,’ 

— + —(sec O1q— tan 41¢ sin A2q 04 /v2) 
2 v1 


, 


b — 
+ ——(sec 4)» — tan 4)» sin 42 v1/v2) 


v1 
which simplifies to equation (11). 

In reducing the results of a travel-time survey, equation (10) is first used in the 
way described previously to yield a first approximation to the time-terms and ve. 
If the results indicate structures steep enough to render the simple theory suspect, 
a second approximation can be made as follows: 


Plot the structure indicated by the first approximation to obtain estimates of 
the angles 62, and 62, for each connection. Also estimate A’ for the various con- 
nections. Finally, we wish to eliminate the angles 4)q and 6)», which will vary from 
one connection to another. This we do by replacing t by t’, where 


AN’ Ha’ cos@ Hy’ cosé 
t' = —+ + ———. (12) 


v2 Vi Vl 
H,' cos @ / cos 644 


We therefore multiply our first estimates of the time-terms by factors of the 
form [(cos@;q/cos@)—1] and use the results to obtain the correction terms t—f’. 
The correction factor is given by the equation 

(cos 6;/cos@—1) = [(1 —r? sin?62)/(1 —r?)}'—1, (13) 


3 


where r = }/v2. Values of the correction term for various values of r and 62 are 
plotted on Figure 2. 

Strictly speaking, the Snell’s Law path shown on Figure 1 can only be drawn 
when the surface of the marker is anticlinal. In fact, successful surveys have been 
conducted over synclinal areas, and it has to be assumed that there is sufficient 
curvature or diffraction of the rays in the marker layer to conduct energy into the 
theoretical “shadow zone’. 

If there is a substantial gradient of velocity in the marker layer the ray paths 
will be curved, and the fact that different rays are refracted upwards from regions 
of different velocity will modify the time-terms. The modification is determined 
entirely by the inclination of the ray path at the point of entry to or exit from the 


marker, and may therefore be treated in the same way as the effect of a dipping 
structure. 


We shall now demonstrate the effectiveness of the corrections by means of a 
worked example. It is assumed that six seismograph stations are distributed over 
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Fic. 2.—Correction terms for various values of v/v2 and 92. 








*Shot?' 





. Shof 3 


Fic. 3.—Hypothetical pattern of shots and stations, laid out over a buried, 
spherical dome. 
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a plane area of the ground, which is underlain by a spherical dome of high-velocity 
material which touches the surface under the central station. Three shot-points 
are laid out on the surface, and the horizontal distance from each shot to each 
station is measured on the plan (Figure 3). The velocities in the upper and lower 
media are taken to be one unit and two units respectively. Ray paths through 
the dome were constructed graphically in accordance with Snell’s law, and the time 
of travel was deduced by summing the times which would be required for a seismic 
wave to cover the segments of the ray paths in each of the two media. The travel- 
time data are listed as Aj; and t;; in Table 1. 


Table 1 


Oi; 


be 


160°7 


wm NNN 
at Ws! 


N 
a 


31°7 


—0°'07 30° 23°8 
—0°63 4 21°5 
+0°23 . 29°7 
+0°04 . 24°3 
+0°45 3° 34°2 
—o'o! 56° 52°6 


aut wn 


—0°02 33° 32°4 
+0°45 35°9 
—0'36 29°7 
—o'1o 48°5 
—0°32 59° 53°3 
+0°34 57° 


I 
2 
3 
4 
5 
6 


First Approximation ‘9 = 2°224 5132 = 1°40 
a1 = 3°19 ‘ 8°73 a3 = 12°74 
6b, = o (arbitrary) 2 0°68 bg = 1°69 


ba = 5°46 = 10°37 be = 14°76 


Second Approximation w 
a.}= 2°97 a 
by —o°30 (arbitrary) b’ 
b’4 5°59 b’ 


'9 = 2°04! 5172 = 0°64 
‘2 = g'Io a’3 = 13°70 
2= 0-40 b’3 = 1°48 
5 = 10°79 b’g = 15°88 


We shall now attempt to recover the form of the dome, assuming that the only 
available information, apart from the travel-time values, is the velocity in the 
upper medium, and the existence of the central outcrop. On entering the values 
of Aj; and tj into the travel-time equations, we find v2 = 2-224 units. The 
residuals for this first approximation are listed in the fifth column of Table 1. 
We determine « by setting 5; = 0, and can then calculate all the remaining time 
terms. These time terms are converted into estimates of the thickness of the upper 
layer by multiplying by v1/+/(1 —v;2/v2?). The first approximation to the structure 





426 P. L. Willmore and A. M. Bancroft 


is the envelope of spheres, each of which has a survey point as centre, and the 
apparent thickness of the upper layer as radius. In our example, the number of 
survey points is not sufficient to determine the three-dimensional structure, but 
we can demonstrate the degree of success of the solution by drawing a section, 
in which the survey points are plotted at the appropriate distances from the centre, 
without regard to their distribution in azimuth (Figure 4). 





Fiest op? oximefon 


nterface 


True inferface 


Fic. 4.—The first approximation to a radial section of the buried structure. 


To obtain data for the second approximation we draw the normals from each 
survey point to the approximate interface, and then draw a vertical line to the 
Earth’s surface from the foot of each normal. The ends of these verticals are the 
accented survey points shown on Figure 3. We now find the slant distances along 
the rays in the lower medium by measuring between the accented survey points, 
and allowing for the difference in length of the verticals at the ends of each ray. 
The results are entered as A’;; in the sixth column of Table 1. We also make a 
rough estimate of 62 for each connection, and derive the corrected travel times 
t’ 43. 

On repeating the calculation with the revised times and distances, we find 
v'2 = 2-041 and derive a new set of time terms and the residuals 5°;;._ On inspect- 
ing the results, it becomes clear that the value of « which makes 5; = o does not 
yield the smoothest fit of the shot-point and station time terms to a common 
envelope. Instead, we choose « to give the smoothest fit, and allow a small negative 
value for 5}. The revised depth circles are shown in Figure 5, and it is clear that 
the solutions are converging quite rapidly towards the true value. 


4." The uncertainties in a time-term survey 


One of the objections which has been raised against the time-term method is 
that one pours all of the data into a large pot and that the conclusions are fished 
out mysteriously after prolonged stirring. Those who are familiar with the pit- 
falls of the conventional interpretations of refraction surveys sometimes wonder 
how the errors of observation are sorted out in the course of calculation. 

In answer to this question, we first consider the degree of control on the estimate 
of the velocity. In this determination all the observational equations have been 
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used, but by subtracting the time terms we have removed those portions of the 
travel-times which are consistently associated with particular points of the network. 
The residuals, 54%, represent the remaining inconsistencies of the observations, 
and we start by assuming that all the uncertainties in 1/v will be reflected in the 


values of 5;,. Working on this assumption, we assign a standard deviation o(1/v) 
to 1/v, where 


(14) 


Second approximahon 
interface 


Fic. 5.—The second approximation to the structure. 


As o%(t) tends to a constant value if the number of available observations of 
given quality is very large, we see that the effect of each new observation in reducing 
the uncertainty in (1/v) is determined by its contribution to =d?;;. Expanding 
d;;, we have 


‘ 


j-l i=1 j=l é=1 

1 I I 
dy = Ay-— > Ay-- > Ay+— > > Ay, (16 
: an : ol = : — | 


j=-m ten 


On the right-hand side of equation (16) the second term is the mean distance 
for all stations which observe the ith shot, the third term is the mean distance of 
all shots observed by the jth station, and the fourth term is the mean distance for 
all the observations. If the stations are all close together, the first two terms will 
almost cancel, and so will the last two. If the shots are all close together, the first 
and third terms will almost cancel, and so will the second and fourth. We see 
from this that the determination of the velocity is impossible without adequate 
dispersal of both shots and recording points. 

The foregoing argument has an important bearing on the practice, which is 
still quite common in crustal studies, of deriving a propagation velocity from the 
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travel-time curve of waves in an underground layer, even when all the waves which 
have been observed have come from a single source. The velocity which is derived 
from such data is the apparent velocity of propagation of a disturbance along the 
surface of the Earth, which will differ from the true velocity in the marker layer if 
the magnitude of the time terms is correlated with distance from the source. Such 
a correlation will arise if the waves are recorded along a single line running across 
a dipping marker layer. Even if observations are made in all directions from the 
source, there will be a correlation between the time terms and distance if the source 
is situated near the centre of a dome or basin in the marker layer. In such cases, 
the trend of the time terms is absorbed into the apparent velocity, and the residuals 
include only experimental errors and the departures of individual time terms from 
the mean trend. The magnitude of any such trend can only be estimated from 
independent structural evidence. In many cases, it is better to assume a plausible 
value for the true velocity in the marker and to discuss the resulting pattern of 
time terms, than it is to use the apparent velocity without further discussion. 

Returning now to the general case, we see from equation (2) that a timing error 
at one station, affecting all shots equally, will be absorbed into the time term of the 
station without showing up in the residuals. In the same way, an error in shot 
timing will be absorbed into the time term for that shot. More subtle errors can 
be introduced by local variations in the velocity of the marker layer, when they are 
distributed so as to affect equally all the travel times which relate to a given shot 
point or seismograph. The time terms for points outside the main area of the 
survey are particularly liable to contain errors of this type, for the seismic rays 
leading to or from any such point w..il all lie close to each other, and will all pass 
through parts of the marker layer which may not be in the path of rays between 
other points of the network. If these parts of the marker layer have an anomalous 
velocity, the anomaly will not receive its proper weight in the velocity equations, 
but will introduce a concealed error into the time term of the outlying survey 
point. A simple example of this principle arises in the case of a line of geophones, 
shot from both ends. Under the geophone spread, where the waves travel in both 
directions, local variations in the velocity of propagation in the marker layer can 
be distinguished from variations in the thickness of the overlying materials, but 
the corresponding separation cannot be made for the materials below the shot 
points. 

It is clearly possible to invent anomalies in the marker layer, which might be 
so distributed in relation to the layout of the observing network as to prevent the 
detection of the errors, but it is much more likely for the concealment of the errors 
to be incomplete. In this case, a critical examination of the results may lead to 
some improvement. The recommended procedure is to write out the residuals 
54; in order of increasing distance Ajj. Negative residuals occurring at the longest 
and shortest distances will suggest that the propagation velocity in the marker is 
increasing with depth. Residuals of a given sign, associated with ray paths of differ- 
ent length, all passing through the same area of the ground, will suggest a velocity 
anomaly of limited horizontal extent. 

The foregoing discussion has revealed an important dilemma in the interpreta- 
tion of refraction studies. If one starts with any given array of sources and detectors 
(including the special arrays, such as reversed profiles, which are used in refraction 
prospecting) one has the choice of determining the velocity and all the time terms 
separately, or of fitting a single travel-time line to all the data and regarding the 
scatter of the time terms as part of the observational error. In the former case one 
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finds, as we have shown above, that certain departures from the assumed type of 
structure can lead to erroneous results, and that the existence of such errors may 
not be indicated in the residuals. In the latter case, one fails to extract all of the 
available information from the primary data, and the unused information increases 
the residuals and leads to an unnecessarily high standard deviation. It appears 
that the way out of the difficulty must be sought in the development of a more 
sophisticated error theory, and the authors would be very glad to hear of any pro- 
gress which can be made in this direction. 

In the light of the foregoing theory, the following working rules are suggested : 


(1) Both the shot-points and receiving locations must form extended patterns, 
to provide adequate determination of the velocity in the marker layer. 

(2) Survey points at which the errors must be closely controlled must be con- 
nected to other survey points in several different directions. Thus, if the time 
terms from a set of seismographs are of vital importance to the survey as a whole, 
the region which contains them should be bracketed by shots. Similarly, vital 
shot points should be bracketed by seismographs. ‘Time terms for the fringe points 
of the survey, where these conditions cannot be satisfied, will be less reliable. 

(3) Do not move stations too frequently between shots, as the maximum weight 
is achieved by observing the greatest possible number of connections within a given 
array of shot-points and seismographs. * 

(4) Arrange for some shot-points to be occupied by seismographs when other 
shots are being fired, so as to enable the arbitrary constant « to be determined. If 
this is impossible, place some of the network points on outcrops of the marker 
layer, or at least at points where the structure down to the marker layer is known. 


5. The historical development of the concept of delay times 


The necessary preliminary to the development of the idea of time terms or 
“delay times’’ is the splitting of the expression for the travel time into two parts: 
the first term denoting the time which the wave would have taken had it travelled 
directly from shot to receiver at the velocity of the high-speed marker layer; and 
the second term expressing the additional time needed to reach and to emerge 
from the high-speed layer. The second term in this expression is familiar as the 
“intercept time” in the case of two plane layers with a boundary parallel to the 
surface. 

Although the assumption of plane stratification is unnecessary to the solution, 
it has been allowed to remain in the great majority of earthquake studies down to 
the present day, and it was the applied geophysicists, through their development 
of arc and profile shooting and of the German “‘laufzeitplan’’, who broke away 
from this concept. In so doing, they acquired a somewhat rigid attachment to 
the use of particular geometric arrays of shot-points and seismometers. In small- 
scale field work, simple patterns are easy to lay out on the ground and are, indeed, 
almost forced on the operators by such instrumental factors as the use of geophone 
cables with evenly-spaced takeouts. In a more general theoretical approach, how- 
ever, we must realize that strict conformity with a prearranged pattern does little 
except simplify the arithmetic, and that the pattern can and should be discarded 
whenever it becomes inconvenient. 

An early method of isolating the time terms for an array of survey points is 
described in the work of the members of the Imperial Geophysical Experimental 

2E 
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Survey (Broughton Edge & Laby 1931). These authors considered the travel 
times of refracted waves, passing between three points, A B and C, arranged in a 
straight line along the surface of the ground. If the seismic waves enter or leave a 
supposedly uniform refractor at points L, M, N and P (Figure 6) then the travel 
times along the various refraction paths are connected by the equation 


Tus+T pn = Tatmp+ Tepnc— Tarpc 





Fic. 6.—Ray paths considered in the ‘‘method of differences.” 


The authors made the simplifying assumption that the velocity contrast be- 
tween the overburden and the marker was so great that the paths in the overburden 
could be considered to be perpendicular to the interface. We now recognize that 
this assumption was unnecessary, and that the I.h.s. of equation (17) is just twice 
the time term of the point B. ; 

This ‘‘method of differences” was refined by J. G. Hagedoorn (1959) under 
the title of the “‘plus-minus method”. He points out that an equation equivalent 
to (17) applies without requiring that the detection point B should be on the sur- 
face of the ground, and therefore defines a system of “‘plus”’ lines in the ground, 
each of which is the locus of points for which the right hand side of equation (17) 
is constant. Hagedoorn also draws a set of ‘‘minus”’ lines, orthogonal to the “‘plus”’ 
lines, each of which is the locus of points for which T4zme—Tppenc— Taxpc is 
constant. The spacing of the “minus” lines along the surface of the ground deter- 
mines the velocity of propagation in the marker layer. It is easy to show that 
Hagedoorn’s ‘‘minus’’ values are a special case of the term cjj in our equation (8), 
and hence the “‘plus-minus” method is the same as the time-term solution for the 
special case of a linear array of detectors and two shot-points. However, the fact 
that the array is linear carries the advantage that velocity anomalies in the marker 
layer can be located exactly, so that the incorporation of such linear arrays into 
a time-term network might sometimes be a valuable method of eliminating local 
uncertainties. 

It appears that the first attempt to survey an extended area was due to Jones 
(1934), whose geometrical arrangement combined the use of arc and profile- 
shooting. For example, shots fired at points S;—S; (Figure 7) were observed on 
the arcs A,—As respectively, and in addition observations were made along the 
line joining the shot-points. These auxiliary shots established the velocity of the 
indirect wave and removed the delay time at the shot-points from the computations. 
In fact, Jones reduced all his times to a common datum (Shot-point No. 1) and then 
plotted a map of the time contours (sum of delay time at each station plus delay 
time at shot-point No. 1). He transformed time contours to depth contours with 
full appreciation that the expressions of the form (H cos 6)/v were characteristic of 
each survey point. 

Jeffreys, meanwhile, independently extended the derivation to the case where 
the depth to the boundary is variable (Jeffreys 1935). He pointed out that although 
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the ray path used in the derivation is not a strict refraction path, since the latter 
is a minimum path, the time of transit will be stationary for small variations 
from it. 

He also showed, by means of an elementary example, that the use of points 
common to both shots and observation renders the interpretation unambiguous. 
In summary he stated that his argument showed “‘that the usual methods are more 


directly related to the depth than might be supposed from the form in which the 
subject is usually presented”’. 


A3 


Fic.7.—Pattern of shots and detectors for combined arc and profile 
method (after Jones). 


As a further modification, Gardner (1939) suggested the use of the “offset 
position” of the survey point, which is the point on the ground, vertically above 
the point at which the seismic ray enters or leaves the marker layer. In drawing 
up plans for areal surveys, Gardner suggested that the survey points should be 
moved for each observation, so that the offset points could be kept constant for 
rays travelling along different azimuths. The present authors feel that this was a 
retrograde step, for we have seen that the calculations required for working with 
fixed survey points are simple and accurate. In a large-scale field operation, it 
would be highly inconvenient to occupy all the survey points needed to hold the 
offset points in constant positions, and there is no guarantee that it would be physi- 
cally possible to do so. Moreover, since materials near the surface of the ground 
are usually more varied than those at depth, the movement of the survey points 
might introduce more errors than it would eliminate. 

In the light of all these developments, it is remarkable that the essential prin- 
ciples of time-term surveys have still not been universally recognized. Unfortu- 
nately, field studies are often carried out without such recognition. As an example 
of the consequences of this fact we shall discuss a paper by one of the present 
authors (Willmore 1949) which describes an operation to which all of the four rules 
could have been applied, but where in fact all four were broken. 

The paper describes the observations which were made on a number of ex- 
plosions in North Germany, the largest of which was at Heligoland, and the 
remainder of which were fired near Soltau. All the explosions were well timed at 
the source, and were observed by mobile field parties. 

The Soltau explosions were fired first, and the resulting waves were observed 
along 4 radial profiles extending about 50km from the source. Preliminary plots 
indicated that waves through the sediments (P;) and through the basement (Py) 
could be distinguished. When the Heligoland explosion was fired, both P, and 
P,, waves were recorded at portable stations at distances up to 314km, and at 
permanent stations as far away as Puy de Déme. 
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In the final discussion, time term theory was applied twice, first to enable 
short-range data to be used in eliminating the effects of overburden, and then in 
an attempt to derive time terms which would describe the passage of P, waves 
through the consolidated sediments. It was at this stage that the shortcomings 
of the operation became evident. If rules (3), and (4), had been obeyed, seismo- 
graphs to record waves from Heligoland would have been set up at the Soltau 
shot-point, and at least some of the points at which the Soltau Py waves had been 
observed. These would have served to determine the velocity of Pj, and an un- 
ambiguous set of time terms. Secondly, a determined effort should have been made 
to record P, from Soltau, for this would have given an unambiguous P,, velocity 
and a further set of time terms. In fact, the networks which recorded Heligoland 
and Soltau waves had no points in common, so that each operation had to be treated 
by itself. Each network, therefore contained only one shot-point in violation of 
rules (1) and (2), and the true velocities in the refracting layers became indetermi- 
nate. This fundamental indeterminancy rendered it impossible to decide whether 
the waves which had been called Pg had all been transmitted through the upper 
part of a single-layer crust, or whether some had been refracted through a basaltic 
layer at greater depth. Finally, the quoted uncertainties in the travel times re- 
presented only the deviations from the times which would have been expected if 
a system of plane, uniform, but possibly dipping refractors had underlain the 
various profiles. As such, the residuals bore no simple relation either to the errors 
of observation, or to uncertainties in the estimates of the true velocities in the 
refracting layers. 


Dominion Observatory, 
Ottawa, 
Canada: 
1960 August 8. 
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The Revision of Earthquake Epicentres, Focal Depths 
and Origin-Times using a High-Speed Computer 


B. A. Bolt 
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Summary 


A program for an automatic computer has been developed to re- 
vise rapidly, provisional foci and origin-times of normal and deep-focus 
earthquakes. For each earthquake, up to 300 equations of condition 
found from P, pP and PKP observations are solved by least squares to 
give a correction to the trial location. Special attention is given to the 
weighting of observations and factors affecting convergence. 

Features of the program are that the theoretical travel-time tables 
are stored in complete form, and after each iteration a list of stations 
with corresponding distances, azimuths, and residuals as well as the 
root-mean-square error is printed. Applications to a 1954 hydrogen 
bomb explosion and a number of earthquakes are described. The re- 
sults suggest that the program may be useful to research organi- 
zations requiring either regular or special location of epicentres. 


1. Introduction 


Location of the foci and origin-times of significant earthquakes is at present 
carried out regularly by several organizations in various countries. In addition to 
this work individual seismologists often re-examine published locations or deter- 
mine locations ab initio from lists of travel-time data; a detailed discussion of the 
main procedures followed in these “special studies” is given by Jeffreys (1959, 
pp. 68-69). Several methods of revision are in common use; the one adopted in 
this work is that of Jeffreys (1959, p. 63) and in essentials is that followed for the 
International Seismological Summary (Hughes 1930, 1933). It involves fitting 
observations of P and PKP travel-times to times in the Jeffreys—Bullen tables 
by least squares. The procedure usually followed in manual computation is to 
reduce the equations of condition to 5 or 6 by combining groups of stations having 
only small variation in distance and azimuth. 

High-speed computing machines have been increasingly applied to data process- 
ing and numerical computation in geophysics (Landsberg 1958; Bullard 1960); 
the aim of the present work was to take advantage of machine speed and storage 
capacity while retaining checks often supplied by individual judgement during 
manual calculation. 

The machine program described, called ID 1, has been tested successfully on 
a wide variety of earthquakes. Three examples outlined in Section 4 show the 


* Received in original form 1960 June 24. 
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revisions of trial locations of a U.S. hydrogen bomb explosion near Bikini in 1954, 


a deep focus earthquake (1934 June 29) and a normal depth earthquake (1957 
March 10). 


2. Formulation of the numerical methods 


In the present work the observational data for the revision of foci are restricted 
to observed arrival times of P, pP and PKP phases. 

Suppose that a preliminary trial has determined a position and origin time for 
an earthquake. Relative to this determination let the true focus be x (deg) to the 
north, y (deg) to the east and 6338zkm deeper; let the true origin time be 7s 
later. Then, for a station at distance A (deg) and azimuth « (deg) from the trial 
epicentre, the equation of condition is, to the first order in x, y and 2, 


an(A)e | 21(8)s - - 


é oz 





T—(x sina+y¥ cos «) 


where ¢ is the theoretical P or PKP travel time to distance A, and € is the P or 
PKP residual (O-C). A sufficient condition for the neglect of terms O(x?, y?) in 
(1) is that |x|+ |y| < [/(2)—1]A. 

Data from n observing stations provide a set of m linear equations of condition 
of type (1). Observational and copying errors as well as other factors require that 
a method of weighting the observations be incorporated before a least square 
solution is made for the corrections 7, x, y, 2. 

Suppose that the true epicentre is to the north-east of the trial one. Then, on 
the average, residuals from stations in the north-east quadrant will be negative 
and those from stations in the south-west quadrant positive, so that continuous 
weighting about the mean residual in at least each quadrant is required. The 
calculation of unrepresentative means, through the presence of chance large 
residuals, is prevented by the rejection of |€| > 40s before the calculation of the 
mean. Generally it follows that such observations receive zero weight. 

The weighting function adopted in ID 1 is that introduced by Jeffreys (1948) 
and found applicable to many distributions of travel-time residuals. Such distribu- 
tions appear to consist of two components, the dominant one approximately 
normally distributed and the other occurring as a background of fairly constant 
density. 

Suppose that the appropriate frequency function is 

in —m)2 
f(é) = Pe. 3 il eee 
o(27)* | 202 | 
where a and s are constants. For a sample of n residuals the maximum likelihood 
estimates for the mean and variance (Jeffreys loc. cit. 1948) are 


= > w(Ei)Ei/ > w(&), 
> w(Ei)(i—m)?/ > w(Ei); 


and 


the weighting function is 


w(f;) = 1/[1 +p exp{(Ei—m)?/20°}], 
where 
fe = o(27)' x as/(1—s). 
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From (5), 4 may be estimated directly from the observations by dividing the 
frequency, as, of large residuals by the difference between as and the frequency 
of residuals near the mode. Strictly, values of ~ and o should be estimated for 
each earthquake but from a number of trials as well as previous numerical work 
(Jeffreys 1936; Stoneley 1939; Bolt 1959), the values x = 0-02, o = /(10) seem 
generally satisfactory; in an automatic iterative method severe weighting near the 
mean is not essential and in the present problem may lead to difficulties. ‘The 
adopted parameters yield the following illustrative weights: 


\fi—m|s 0 2 4 6 8 10 12 
wi) o98 o9¢8 o96 ogo O67 025 0°04 


If the appropriate diagonal weighting matrix W is formed from 4/(w(£;)|, the 
n weighted equations of condition derived from equations of type (1) may be 
written in matrix form as 


WAx = WB (6) 


where A is the m x 4 matrix of coefficients of 7, x, y, z and B the x 1 matrix formed 
from the residuals £;. Equation (6) is stored and manipulated within the machine 
in the form 


Cy = 0, (7) 


where C = (WA|WB) is the augmented mx 5 matrix, and y is the vector ( *,). 
The normal equations are then 


Dy = 0, (8) 
where 


D = ATW'C. 


The epicentral distances and azimuths of the observing stations from the epi- 
centre are calculated by the formulae introduced by Turner (1915), after computa- 
tion of the geocentric direction cosines from the co-ordinates. Ellipticity adjust- 
ments are computed using the approximate formula derived by Bullen (1937). 
Elements of the first four columns of C are generated by short subroutines. From 
stored theoretical travel-time tables, the derivatives ¢t/0A and 0t/é@z in (1) are 
formed from the first differences, and the residuals by two-way interpolation. 

Solution of equations (8) is by Gaussian reduction combined with searching 
columns for the maximum element before pivoting. All matrix elements are 
reduced to the same order of magnitude by dividing the elements Cj4 in (7) by 
100 and scaling the resulting corrections to focal depth accordingly. 

In the case of normal earthquakes, the fourth element of the vector x is put 
equal to zero and only the first column of the travel-time tables is entered. Deep- 
focus earthquakes are distinguished by the presence of both P and pP travel-times. 
Initially pP-P intervals alone are used to estimate a mean focal depth and hence 
the column of the travel-time table to be entered, but in the least-squares calcula- 
tions the pP-P data augment the P and PKP data by the addition of equations 
of type (1), modified by placing the first three terms zero. 

The solution provides values of 7, x, y, z which are used to revise the trial 
origin-time, epicentre and focal-depth and the procedure is repeated. After each 
iteration the root-mean-square or standard error (se) relative to the new weighting 
function w(&), given by 


se = [> w(Eié?/> w(&)), (9) 
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is computed for the residuals. As the four mean residuals, m, approach zero, (9) 
approximates to the standard deviation defined by (3), in agreement with the 
assumption that the true focus and origin-time are such that the residuals for each 
of the stations could be reduced to zero except for the presence of independent 
random errors in the observations. 

In practice there is no assurance that the process described converges to the 
true solution for a particular earthquake. Iteration may fail, for example, if observa- 
tions are badly inconsistent among themselves (as in multiple shocks). For single 
shocks reliable locations depend upon the degree of fulfilment of two criteria: 


(1) observations should be evenly distributed by azimuth; and 
(2) within each quadrant, stations should be evenly distributed by distance. 


Jeffreys (1959, p. 64) estimates that standard errors in the corrected latitude and 
longitude of about 0°-o5 are attainable if the earthquake is well-observed in this 
sense. A useful measure of convergence is the trend of the magnitude of the 
standard error but its actual value may underestimate the uncertainty in the loca- 
tion if there is correlation between the errors at groups of stations. Convergence 
of the method in practice may be inferred from examples discussed in Section 4. 


3. Description of the program ID 1 


ID 1 has been written for an IBM 704 data processing machine which executes 
instructions at the rate of about 40000 per second. A maximum of 300 P and 
PKP observations can be used for the revision of each earthquake and any number 
of earthquakes can be revised in a single production run. The program, travel- 
time tables and observational data require 16000 memory locations. The seismic- 
data input which is on punched cards is of three kinds: 

Travel-time data. ‘To facilitate linear interpolation the standard Jeffreys— 
Bullen tables for the main seismic phases give travel-times for 14 depths, namely, 
at the surface and at 0-00 R(o-o1 R)o-12 R, where R = 6338kms. In a machine 
with spacious memory it is convenient to store the times in this form rather than 
generate them from functions previously fitted to tabular values in various distance 
ranges. Consequently the P travel-time tables (0° < A < 105°) are stored as an 
106 x 14 array, the pP-P tables (25° < A < 100°) as a 76x14 array; and the 
PKIKP tables (120° < A < 140°) as a 21 x 14 array. Tables for ellipticity adjust- 
ments are also stored in the memory. 

Earthquake data. 'The date, trial co-ordinates, and origin-time are punched 
on a single card; no trial focal depth is punched. The format of this card defines 
the printing arrangement of the revised focus and origin-time (cf. Section 4). 

Station and observational data. A library of cards has been prepared containing 
all seismological stations listed by the J.S.S. (1951) and by the Bureau Central 
International Seismologie. The cards are punched as follows: 


RIVERVIEW LAT —33:83 LONG+ 151-15 P PP PKP. 
After punching of travel-times on copies of relevant cards, they are assembled in 
the order required for printing of the station names. 
A flow diagram on the opposite page summarizes the structure of the program. 
Iterations stop when |x|, |y| < 0°-1 and || < o-1s, or when the iteration counter 


equals a specified number NA; usually two iterations are sufficient to reduce the 
se to about 2s. 





Earthquake epicentres, focal depths and origin-times using a high-speed computer 437 


Flow Diagram 


Read into computer Jeffreys—Bullen tables 
pa Read in the assumed position and origin-time of shock 


Read in iteration counter = NA (NA = 3 suitable) 


Read in additional station counter = NB 


v 
p——PRead in station names, positions and travel-time data 


\—2 Compute distances, azimuths and ellipticity corrections 
(after first passage) 


<< 


«< Assigned transfer 


| (first passage only) 
Test for focal depth (presence of pP times) 





If shallow, jump Sub-routine II 


Estimate focal depth from pP-—P intervals 





LepCalculate P, PKP residuals 


Y 
Find appropriate weights treating each quadrant separ- 
ately 


v 


Calculate the root-mean-square error 


v 
Print out station names, distances, azimuths and residuals 


° 





Test if NA 
4 ° 
Test if NB ° 
° 











N Halt or proceed to next shock 





Form weighted equations of condition and solve by least 
squares 


Correct latitude, longitude, focal depth and origin-time 
of shock 


v 
Revise relevant observed travel-time data 








a 
= 


y 
Test magnitude of corrections 





Input I 


Input IT 


Sub-routine | 


Sub-routine I] 


Sub-routine III 


Output 


If jump, put 

NA = NA=1 

If jump, put 

NB =o 
(additional stations 
read in) 


Stop 


Sub-routine IV 


Continue or stop 
iterations 
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Station names, distances, azimuths and residuals relative to the current solu- 
tion are printed out after each revision and after the final iteration this information 
is printed also for a pre-specified number of stations, additional to those used in 
the revision. Certain other information can be printed at the end of each iteration 
by the use of sense switches. This includes a list of geocentric direction cosines 
and the elements of the matrices C and D. A feature of the output is that results 
are listed in similar form to that in the J..S.S. and are suitable for printing without 
further transcription. 

Revisions for over 30 earthquakes, each requiring an average of three iterations, 
indicate a computation time under one minute per earthquake. 


4- Application of ID 1 to a nuclear explosion and two earthquakes 

An early test computation used 19 P observations (Kogan 1960) from the 
nuclear explosion on 1954 February 28, in the vicinity of Bikini atoll. The trial 
and first-revised positions, origin-times (to), and standard errors were: 


Trial 11°-0O0N 165°-00E, 18h 45m §s, se = 8-94s 
Iteration I 11°°73N 165°-32E, 18h 44m 57-798, se = 1°258. 


The printed output after the second (and final) iteration is reproduced below: 


H BOMB BIK FEB 28 1954 1le74 165632 FD=-3340 18 eH 440M 572805 
SE= 1425 


1D1 


DELTA(DEG) AZ(DEG) RESISEC) 
NOUMEA 33285 176213 0.68 
YUZNOSAKL 39294 335279 2048 
BRISBANE 40.77 196284 049 
VLADIVSTK 42259 323442 1.69 
BAGUIO 43.58 281249 -0262 
RIVERVIEW 47228 196.611 -0.09 
KIAKHTA 61.35 321234 -1206 
IRKUTSK 63218 322-97 -0225 
VICTORIA 68.88 42245 -1.53 
PASADENA 72239 57279 0.93 
TINEMAHA 7223 54270 2200 
SEMIPALAT 77293 319-230 ~jo32 
ALMATA 80.89 31231 -1.250 
FRUNSE 82266 312218 0266 
TASHKENT 86.83 311239 0.88 
SAMARKAND 88.92 310220 0-09 
QUETTA 91-21 301-03 -0.70 
FAYETTVIL 92.00 54.73 -2034 
UPSALA 104,08 343.485 0247 


(Since the J.-B. tables conventionally take the base of the crust as 0-00 R, shallow 
focus sources are printed with FD = —33:okm.) The officially stated position 
and origin-time are 11°-69 N 165°:27E, to = 18h 45m 0-00s. Using 232 P obser- 
vations with 33° < A < 104° from 9 U.S. nuclear explosions in the Pacific, Kogan (loc. 
cit. 1960) gives the mean deviation from the J.-B. tables as —1-8+0°6s. Applying 
this as a correction to fo, the official time is duplicated within 0-5 s, a result in agree- 
ment with a hand computation by Mme Labrouste (1958) who finds to = 18h 
44m 58-2+0°5s from slightly different data. 
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The Celebes Sea earthquake of 1934 June 29 allows a test of ID 1 in the case of 
very deep focus. Equations of condition from 55 P observations for 5° < A < 98 


2 PKP observations for 123° < A < 136°, and 13 observations of pP adjusted a 
trial epicentre as follows: 


Trial 6°-10S 123°-40E FD = 668-1km to = 8h 25m 20-008 se = 4°658 
Revised 6°-91S 123°-23E FD = 649-5 km to = 8h 25m 13-778 se = 2°448. 


The trial epicentre estimated originally by Berlage (1936) and adopted by the 
I.S.S., has been revised previously by Jeffreys (1942). He combined the normal 
equations from P, S and SKS and obtained a similar solution, namely: 


6°-91S + 0°-03 123°-19g E+ 0°04, FD = 647-2km+ 4-4, 
to = 8h 25m 13°08s+0-4 


As a preliminary to a study of travel-times of core phases from 10 Aleutian 
earthquakes selected from the series in 1957, the determinations listed in the 
B.C.I.S. were re-examined using ID 1. The successive corrections to the trial 
epicentre of the least well-determined shock, 1957 March 10, show quantitatively 
the convergence in such a case: 


Trial 52°-00N 176°-00W to = 3h6m 2-00s se = 10°78 
Iteration I 51°-60N 174°-78W 3h6m 987s 3°58 
Il 51°-72N_ 174°:28W 3h 6m 10°36s 2°18 
III 51°-76N 174°10W 3h 6m 10°58s 1°gs 


After the third iteration, except for residuals of 27s (Seattle), —19s (Apia) and 


11s (Tacubaya), the residuals for the 37 stations used are within +38, in agree- 
ment with the standard error. 


5. Other programs 


An extension to ID 1, designed for the revision of travel-time curves, has also 
been programmed for the IBM 704. This program requires 22000 memory posi- 
tions. Suppose m is the number of earthquakes used. After the final iteration for 
each earthquake, the revised focus and origin-time are used to calculate, for the 
phase under study, weights and adjusted travel-times to as many as 150 stations. 
These are stored and accumulated with similar values from each of the remaining 
m—1 earthquakes. Polynomials of degree up to a predetermined value are then 
fitted to the weighted travel-times and root-mean-square errors, indicating the 
goodness of fit, calculated in each case. 

Three other programs for the automatic location or revision of provisional 
foci are known to the author. At the 1960 Vancouver Meeting of the Seismological 
Society of America, Brazee and Gunst re: -ted progress using Geiger’s method 
(Geiger 1910) with an IBM 650 computer at the United States Coast and Geodetic 
Survey, Washington*; a similar method has been coded for SILLIAC at the 
University of Sydney, Australia, by Dr Leo Howard. Finally, Dr Eugene Herrin 
at Dallas Seismological Observatory, Texas, kindly outlined in correspondence a 
program, similar in many respects to ID 1, coded for UNIVAC 1103. Experience 
with ID 1 indicates that such programs may be of considerable value to organiza- 
tions regularly determining earthquake epicentres, as well as to individual research 
workers requiring as precise determinations of earthquake elements as possible. 
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Program ID 1 was developed on the IBM 704 computer at the IBM Research 
Centre at Poughkeepsie, New York, The author carried out the work under U.S. 
National Science Foundation Contract G-3485 while a recipient of a Fulbright 
Research Scholarship Travel-Grant. 


Lamont Geological Observatory, University of Sydney, 
(Columbia University), Australia. 
Palisades, New York. 
1960 June 20. 


(Lamont Geological Observatory. Contribution No. 442). 


* Note added in proof.—Since the card for 1960 August 8, the U.S.C.G.S. 
preliminary epicentres have been determined using an electronic computer. 
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Research Note 


Experimental Investigation on the Properties of 
Stoneley Waves 


R. J. Donato 


(Received 1960 June 9) 


Method 


When a solid has a plane free boundary, a Rayleigh wave can be propagated 
along the surface. If the solid is covered with a liquid the Rayleigh wave still 
exists, but is attenuated, and another type of surface wave, the Stoneley wave, 
appears. Its velocity is fixed by the clastic properties of the solid and the liquid, 
and its amplitude falls off exponentially in a direction normal to the boundary 
both in the liquid and the solid. This note describes measurements made both on 
the velocity and on the attenuation of the Stoneley wave and compares the results 
with theoretically predicted values. 

A barium titanate transducer placed on a paraffin wax block immersed in water 
excited the Stoneley wave. A similar transducer, whose height above the wax 
block could be varied, received the wave. Both transducers were backed to reduce 
“ringing”. ‘The received signal was displayed on an oscilloscope together with a 
time “ruler’’ marked at 1oyusec and r1oopsec intervals. An arrangement was 
included to delay the trigger to the oscilloscope so that the oscillograms shown in 
Figure 1 do not start from zero time. These oscillograms show the different times 
of arrival of (a) the Stoneley wave when the block is covered with water, and (5) the 
Rayleigh wave when the block has a free surface and both transducers rest on its 
surface. 

Time-distance curves were plotted for the dilatational and Stoneley waves for 
the wax block immersed in water. From this information we determine the veloci- 
ties of these waves. With the paraffin wax block in air we can determine the 
Rayleigh wave velocity. The dilatational and Rayleigh wave velocities are sufficient 
to specify the elastic constants of the paraffin wax, and with the added values of 
the velocity of sound in water and the densities of paraffin wax and water we can 
calculate the Stoneley wave velocity. 


Results 


The velocities of the dilatational and Stoneley waves were measured as 
21°3x 104cm/s and 6-5 x 104cm/s, and the Rayleigh wave 8-3 x 104cm/s. If « 
and £ are the dilatational and shear wave velocities in paraffin wax, ao the 
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velocity of sound in water, po, p; the densities of water and paraffin wax, then the 
velocity C of the Stoneley wave is given by 


(2B®/C®— 1)? — 4(8%/C2) ¥/(B2/C2 — B2/a2)y/(B2/C2— 1) + 
+ (po/pi) V/(B?/C2— B2/22)//(B2/C2— BP/ag2) = 0 


(see, for instance, Ewing, Jardetzky & Press 1957). 


(1) 


ee 


3008 400s 500ms 


(a) 


Fic. 1.—Oscillograms of (a) Stoneley, and (6) Rayleigh waves. 


Knowing the Rayleigh velocity and the dilatational velocity we find o to be 
0°40 from tables (Knopoff 1952). ‘This fixes B and so we can solve equation (1) 
with % = 14°8x 104cm/s and po/p; = 1°08. The solution of the equation was 
performed numerically and gave the Stoneley wave velocity as 6-8 x 104cm/s 
compared with the measured value of 6-5 x 104cm/s. 

The theoretical attenuation factor in the water normal to the interface is 
ky/(1 —C?/ao”) nepers/cm. The wave number k is determined from an estimate 
of predominant frequency in the Stoneley pulse shown in Figure 1. Substituting 
for k, C, ao, we get the attenuation factor to be 3-00 nepers/cm. Figure 2 gives 
the results for two sets of determinations on the amplitude of the first peak of the 
Stoneley pulse. A least squares fit gives an attenuation constant of (3:32 +0-06) 
nepers/cm. This is significantly different from the calculated value of 3-00; the 
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difference chiefly arises in the estimate of the frequency to be used in the theoreti- 
cal expression, as the attenuation is directly proportional to this frequency. More- 
over, the spectrum of the pulse, and consequently its predominant frequency, 
varies with different heights above the interface. 
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Fic. 2.—Amplitude variation of Stoneley wave with height above interface. 
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Report on Progress in Geophysics 


Palaeomagnetic Directions and Pole Positions, Part II 
Pole numbers 2/1 to 2/41 and 1/71 (m1) 


Compiled by E. Irving 
(Received 1960 July 11) 


This list incorporates the results available to the author up to April 1960. Cases 
are now arising where old results are being modified. Changes may arise because 
of the reworking of old primary data or by the incorporation of some new data from 
later sampling. Such results retain the old number followed by (m 1), (m 2), etc., 
signifying the first, second such modification, the nature of the change being 
indicated in footnotes. Where a rock formation which has been previously studied 
is re-sampled so that the primary data are entirely new, the result is given a new 
list number. 

The symbol given in column 3 refers to the geological age of the rock unit studied 
and not to the directions of magnetization and pole positions unless evidence (noted 
in column g) is described in the original which allow it to be supposed that the direc- 
tions of magnetization were acquired at the time the rock unit was formed. It is 
hoped that more evidence of this type will be forthcoming in the future. Results 
calculated from rock units which do not cover in time the full range of the secular 
variation will be subject to errors which are not always taken account of in the pub- 
lished errors. At the moment there is no general agreement about the method for 
estimating errors, so that although the values listed in column 15 are given at the 
same probability (P = 0-05) they are not always strictly comparable as different 
schemes for weighting specimens, samples and collecting sites have been used. 
For details of these readers are referred to the originals. These last remarks are 
made in order that a not too uncritical use of the list data will be made. 

The results are compiled without prejudice and are intended only as a list of 
what has been done in palaeomagnetis:m and as a guide to the literature. The 
inclusion of a result in this list does not necesarily mean that the compiler agrees 
with the methods used. Because of the time required to abstract data, which are 
often presented in widely differing forms in periodicals in many languages, it will 
not always be possible to include results in the list immediately following their 
publication. However, all efforts will be made to ensure that these lists are up to 
date, and to help in this regard it would be very greatly appreciated if notification 
of new work, comments and criticisms could be sent directly to the compiler. 


Dept. of Geophysics, 
Australian National University, 
Canberra, A.C.T.., 
Australia: 


1960 July. 
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Sampling 


| 
| 


Country No. Age Rock unit 


Thickness in m 


r 


(x Q Baku stage, sediments, Cheleken (39-5 N, 53 E) 
Q Khazar beds, sediments, Cheleken (39°5 N, 53 E) 


[ Cheteken (39°5 N, 53E) ‘= zone 


N 
Po) 
ce) 


Ir. zone 
3 Tp-Qp* Apsheronsk 
sediments 1 Kiuendag (39 N, 55 E) up. zone 
Ir. zone 
Cheleken (39-5 N, 53 E) 
Small Balkhan (39-5 N, eso 


reversed 


4 Tp-Qp* Akchagylsk 
sediments 4 55 E) 
| Kiuendag (39 N, 55 E) 
U.S.S.R. Mean result for Apsheronsk and Akchagylsk stages 
(Turkmenia) Cup. zone 
Cheleken » third zone 
Tp Red-coloured ( (39°5 N, 53E) °* second zone 
sequence, Llr. zone 
sediments Small Balkhan faundl 
(390°5N, 55E) reversed 


Small Balkhan normal 
(39°5 N, 55 E) reversed 
;reat Balkhan normal 
L L (39°5 N, 54°5 E) reversed 


6 Tpa-o Sediments ) c 





{ 7 Tp(u) _ Basalts, Tannu Basin (35 N, 137E), two flows 
, 8 Ter’ Andesites, Mitaki, 3 sites, 2 flows 

apan (Honshu post ‘ : fe ; : 

Japan memnaned 4 g ‘Tm(m-u) Two flows (andesite) and 2 sheets (dolerite and andesite) 
(ro Tm () Two flows (andesite) and 1 dolerite sheet 


Antarctica 1 T-Q Volcanics, Hallet Cove (72S, 171 E) 


England 12 To-m_ Basic dykes, 7 sites, 4 dykes (55-5 N, 3 W) 200 km 


Scotland Te-m < Gabbro, Ardnamurchan, 3 localities (56-7 N, 6:2 W) ~5km 


| Layered gabbro, Rum (57-0 N, 6-4 W) 500 


14 
L415 


{x} ( Lavas and intrusives, Skye, 53 units (57-4 N, 6-3 W) ~ 1000 4°0km 


16 | Dykes and sills, 16 units 
Ireland 17> Te? Antrim igneous suite Middle Basalts, 6 flows \ 100 km< 
18 } (55:1 N, 6-4 W) Lower Basalts, 19 flows 3 


Japan 19 Ki-m Inkstone Series, red shls., 3 sites (35°5 N, 138-5 E) 


(20 Ku Sediments, Small and Great Balkhan (39-5 N, 55 E) 
J 21 Kl Lr. Albian sediments, Great Balkhan ese 
U.S.S.R. : (39°5 N, 54°5E) \ reversed 
(Turkmenia) }2. Jm-u Sediments, Great Balkhan (39-5 N, 54°5 E) 
J-K Mean of J and K 
England E- Ji Cotswold Sands, reversed only (51-5 N, 2-5 E), 2 sites 
24 Jl Midford Sands (51-5 N, 2°5 E), 3 sites 


France 25 Ji Lr. Hettangian tuffs, 5 sites (43-0 N, 1-6 E) 





Sampling 
EE 


No. Samples 
No. Specimens 


Ni 


nN 
° 


Stability 


‘Treatment 
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Polarity 
—— 


Mean direction of 


magnetization 


D 


t of 


+40 
+40 


2626 —64:'1 
103°5 +69°8 


549 +58-9 


Pole position 


Long. 


168 E 
160 W 


108 W 


57W 
109 W 
152W 
113 W 


86 W 
76W 
159° 6E 
165°‘0E 
150°9 E 
158E 
109 E 
1i58W 
153W 
164E 
169 E 


165E 
167E 


59°3E 
41°0E 


76:6E 


Reference 
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Country 


Antarctica 


Rhodesia 


England 


Germany 


U.S.S.R. 


Scotland 


Norway 


Antarctica 


E. Irving 


Rock unit 


Lsteral spread 


Thickness in m 


r 


Sills and dykes, 7 intrusions 
Ferrar Ferrar Glacier, 5 sites (78S, 161 E) 30km 
Dolerites \ Wright Valley, 8 sites (77-5 S, 161°5 E) 40km 
Beacon Ss., 3 sites (78S, 161 E) 
10 basic dykes in basement, Wright Valley 
(77°5S, 161-5 E) 


Cave Ss.,*(8) 
Sandstones 


Vitloosian sediments, R. Vyatka (59 N, 50E) 
Tartarian sediments, 3 localities, R. Vyatka 

(59 N, 50°5 E) 
Kazanian sediments, 2 localities, Perm (57-5 N, 55 E) 100 km 
Ufimian sediments, 2 localities, Perm (57-5 N, 56 E) 100 km 
Mean for Up. Permian of Russia 


20km 


Whin Sill, dolerite, 34 sites (55 N, 2 W) 120km 


Four igneous units, sills and lavas, Nahe trough (50 N, 8 E) 40 
Tremadocian sediments (60 N, 30 E) 


Younger Gabbros, Aberdeenshire, 16 sites 2 000 km? 
57°5 N, 2:5 W) 


Sparagmites of L. Mjosen (60-5 N, 11 E) 


Gneiss Ongul Islands, 3 sites (69S, 40 E) 





Palaeomagnetic directions and pole positions, Part II 


Sampling 
Mean direction of 
Polarity magnetization Pole position 
tesa aati 
_ D Long. Error 
dy, dx 


o. Samples 


No. Specimens 
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Reference 
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173E 
179W 
178 E 12, 19 


168-9 E 2°8, 2:8 
174E 
169E 


176E 


nt 7 q 179E 
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Pole calculated by the compiler from data given in the originals cited. Poles calculated in entries 16, 17, 
18 and 32 are in excellent agreement (better than 1°) with those obtained by Cox & Doell (1960) 
Directions corrected for instability where necessary by a process of rotating directions along great circles 
away from the present Earth’s field direction (see Khramov, 1958, p. 93). 
Mean directions obtained as an average of samples. 
The Apsheronsk beds overlie the Akchagylsk but it is a matter of debate as to whether they are Pliocene or 
or Pleistocene. 
Dykes of pre-Beacon Ss. age. 
Directions corrected to gneissic banding. 
Sediments may have been heated and remagnetized at the time the adjacent Ferrar Dolerites were 
intruded. 

(8) Results from the same formation as 1/71 with some new data added 








Report on Progress in Geophysics 


Geophysics in Australia 


J. C. Jaeger and R. F. Thyer 


Introduction 


This review is restricted to the geophysics of the solid Earth and to the main- 
land of Australia, though occasional reference is made to associated regions. 

The Commonwealth Bureau of Mineral Resources, Geology and Geophysics, 
is responsible for gravity and magnetic mapping of Australia and carries out 
extensive work on field and station seismology and other aspects of fundamental 
geophysics: this work has been extended to the Antarctic with the collaboration 
of the Antarctic Division of the Commonwealth Department of External Affairs. 
Geophysical prospecting carried out by other Government agencies and by private 
companies provides further fundamental information for studies of the crust. 

In the Universities, the Australian National University possesses a Department 
of Geophysics; the Universities of Sydney, Melbourne and Tasmania provide 
optional courses in geophysics in their geology departments; the Physics Depart- 
ments of the Universities of Western Australia, South Australia and New South 
Wales, and the Geology Department of the University of Queensland, are all 
active in some branches of the subject. 

Magnetic observatories were established by Ross in Hobart in 1840 and by 
Neumayer in Melbourne over a century ago (Dooley 1958). Seismological stations 
equipped with Milne-Shaw instruments were established at the astronomical 
observatories in Perth and Melbourne in 1go1, and in Adelaide and Sydney in 
1906. Riverview, Australia’s best known seismological station, began operating 
in 1909. 


Gravity 


During the last few years there has been a substantial increase in the numbe1 
of gravity observations made in Australia. More than 10000 new stations have been 
established during the past 12 months, about 70 per cent by private companies and 
the remainder by Government agencies and Universities. 

The Commonwealth Bureau of Mineral Resources, Geology and Geophysics, 
is the national authority for collecting gravity data for regional and geodetic pur- 
poses. The private companies have generously made data available which have been 
used by the Bureau in preparing a Bouguer anomaly map of Australia. A pre- 
liminary version based on average values per one degree squares is shown as 
Figure 1. 

The datum for this work is the National Gravity Base Station (NGBS) at the 
Bureau’s Geophysical Laboratory in Footscray, a suburb of Melbourne, Victoria. 
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International connections to NGBS have been made by various overseas and Aus- 
tralian observers listed in Table 1. Groups of observations by individual authori- 
ties have been adjusted to datum through a network of 59 pendulum stations 
established by the Bureau during 1950 and 1951 with the Cambridge Pendulums 
(Dooley & others 1959). The estimated standard error of the Cambridge Pendu- 
lum measurements relative to NGBS is +0-6 mgal. 

Underwater gravity observations have been made by the Bureau with a North 
American underwater gravity meter in Port Phillip Bay, Victoria; along the north 
coast of Australia from Cape York to Wyndham, W.A., via the Gulf of Carpentaria 
and Darwin; in the Gulf of Papua between Port Moresby and Bramble Cay. A 
series of traverses has also been made across the Great Barrier Reef in Queensland 
from the coast for distances up to 80 miles. 

In 1956, the Bureau co-operated with Columbia University and the Royal 
Australian Navy in carrying out a submarine gravity survey of the south-west 
Pacific Ocean. Vening Meinesz pendulum equipment was used in HM Submarine 
Telemachus along the route Sydney—Wellington—Auckland—Tonga—Fiji—Norfolk 
Island—Sydney, including several traverses across the Kermadec~Tonga Trench. 

Only a few isostatic reductions have yet been made. ‘The Cambridge Pendulum 
stations of the 1950-51 survey have been reduced for Hayford theory with thick- 
nesses of 56-9km, 80km, 96km, 113-7 km; and for Airy theory with thicknesses of 
20km, 30km, 40km and 60km (Dooley & others 1960). A few stations have been 
reduced for Hayford 113-7km only in the Carnarvon Basin, W.A., and in the 
vicinity of Cairns, Queensland. 

No measurements of Earth tides have yet been made in Australia but equip- 
ment comprising a North American underwater gravity meter and recorder is at 
present undergoing tests in the Bureau’s laboratory prior to commencing a pro- 
gramme of observation at Melbourne. Measurements will later be made at the 
Bureau’s geophysical observatories in Mundaring, W.A., and Port Moresby. 

The provisional Bouguer anomaly map (Figure 1) will no doubt be modified as 
more gravity data become available; however, there are a few major features already 
well established and these can be correlated with major geological features. The 
gravity low on the south-western coast near Perth, W.A., is due to a rift-like 
feature containing more than 20000 feet of sediments (Thyer & Everingham 
1956). The major gravity low west of Alice Springs, N.T., can be correlated with 
a great thickness of sediments in the Amadeus Basin (Marshall & Narain 1954). 
It is of interest to note that the Broken Hill, N.S.W., Mt. Isa-Cloncurry and the 
Darwin—Katherine mineral fields are associated with regional gravity highs, pos- 
sibly due to the extensive development of dense Archean rocks. 

There is every reason to expect that the present scale of activities in gravity 
field work will be maintained for at least some few years. Geophysical prospecting 
for oil is more active now than in any previous period. Helicopters are being 
used on an increasing scale for transport; during 1959 an area of 25000 square 
miles was covered on a 10 x § mile grid and plans are under way to survey about 
100 000 square miles with a similar station density in 1960. 


Geomagnetism 


The measurement of the Earth’s magnetic field on a national basis is the 
responsibility of the Commonwealth Bureau of Mineral Resources, Geology and 
Geophysics. For this purpose, the Bureau operates magnetic observatories at 
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Toolangi (Victoria) Port Moresby (Papua) and Gnangara (near Perth, W.A.); 
the last is part of the Mundaring Geophysical Observatory. The Watheroo Mag- 
netic Observatory, which was established over 40 years ago by the Carnegie Institu- 
tion of Washington, was operated by the Bureau from 1947 to 1959. The equip- 
ment and personnel were transferred to Mundaring, about 100 miles further south, 
and the Watheroco Magnetic Observatory was closed down in 1959. Magnetic 
observatories are also operated by the Bureau at the Australian Antarctic Research 
Stations at Macquarie Island, Mawson and Wilkes. 

Hourly mean values from the observatories are published periodically. Those 
from Toolangi for the past decade are being published (van der Waal & Brooks 
1959, van der Waal 1960, van der Waal & Sorenson 1960). Watheroo results have 
not been published since 1947, but it is hoped to bring these up to date soon. 
More current data, such as mean monthly values, K-indices, etc., are published in 
the “Geophysical Observatory Reports”, which are exchanged with institutions 
throughout the world. 

A total of over 700 field magnetic stations have been occupied in Australia and 
New Guinea since 1910. About twenty of these are reoccupied every 5 years to 
determine secular variation. These form the sources of information for isogonic 
maps which are issued every five years; the latest applies to the epoch 1960-5 
(Parkinson 1959a). Isomagnetic maps of all elements for the epoch 1957-5 have 
been prepared and are being printed. 

As well as land stations, the magnetic field has been measured recently over 
the oceans bordering Australia by the Russian survey ship Zarya, the American 
survey ship Vema and the U.S. Navy’s worldwide aeromagnetic survey. 

Large areas of Australia have been surveyed by airborne magnetometer. 
Figure 2 shows the area surveyed to the end of 1959. Although this work was 
designed primarily to help in the search for minerals, a knowledge of the mag- 
netic field in detail over large areas is often of interest to fundamental geophysics. 
The most noteworthy of the aeromagnetic surveys being carried out by the Bureau 
covers an area of over goooo square miles along the west coast of Australia. The 
area contains an interesting sedimentary basin and its contact with shield rocks 
to the east. 

Some investigations of a more specialized nature are being carried out by the 
Bureau in the field of geomagnetism. The following may be of interest: 

(i) The aeromagnetic survey over the deep sedimentary basin near the west 
coast demonstrated that a large area was unusually quiet magnetically. Taking 
advantage of this, an attempt is being made to detect the ‘‘non-potential” compo- 
nent of the magnetic field by direct measurement around a circuit enclosing 
1 000 square miles. 

(ii) About four years ago a correlation between vertical and horizontal variations 
was discovered at Watheroo. It is not yet clear whether this should be considered 
an anomaly, or whether this phenomenon is common to most places in the world. 
Temporary recording sites are being established (using Askania variographs) to 
find the regional distribution of this correlation, which may be connected with 
conductivity discontinuities deep underground (Parkinson 1959b). 


Seismology 

The Australian mainland is an area of relatively low seismicity. Up to 1951 
only 80 epicentres of earthquakes of magnitudes 3 to 7 had been located (Burke- 
Gaffney 1951). ‘These shocks, and a fairly large number of much smaller ones which 
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can be recorded with instruments of high magnification, are ascribed mainly to the 
adjustment of local geological structures (Jaeger & Browne 1958). Because of the 
paucity of local shocks, seismologists have been interested mainly in teleseismology, 
though two earthquakes in the Adelaide region have been studied in detail (Bolt 
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Fic. 2.—Australia and New Guinea, showing airborne magnetic surveys 
to May 1960. 


1956; Bvl'en & Bolt 1956) and several observatories are now studying small local 
tremors and their geological interpretation. 

Of the major seismological stations in Australia, Mundaring and Melbourne 
are operated by the Bureau of Mineral Resources which is also responsible for 
Macquarie Island, and Mawson and Wilkes in Antarctica, as well as Port Moresby 
and Rabaul in the active New Guinea area. The Perth station is at the Observa- 
tory, and Riverview in Sydney is operated by the Society of Jesus. Adelaide, 
Hobart and Canberra are run by their respective Universities. The University 
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of Queensland operates Brisbane and Charters Towers, and in addition has 
specialized in the study of microseisms (Jones 1951) and has stations at Brisbane 
and Townsville for this work: additional stations for this purpose will be estab- 
lished at Mackay and Gladstone by the Commonwealth Bureau of Meteorology. 

Two networks, each of four stations with Benioff instruments, have been 
established by the Snowy Mountains Hydro-electric Authority and the Sydney 
Metropolitan Water Board in the Snowy Mountains and Sydney regions, respec- 
tively, to locate small tremors for engineering purposes. The records of these 
stations are handled with those of the Canberra station and provide an excellent 
coverage of S.E. Australia. Epicentres are calculated on an I.B.M. computer: 
about 15 local tremors being located annually, as well as many smaller ones with 
less accuracy. In Tasmania there is a network of four stations fully telemetered 
to Hobart where recording will be on magnetic tape. The three stations remote 
from Hobart form an equilateral triangle with sides 70 miles in length. At present 
conventional seismometers are being used, but alternating current seismometers 
have been developed (Newstead & Watt 1960) and are to be installed. 

One interesting result of the improvements in instrumentation in Australia in 
the last two years is that the seismicity of the mid-oceanic ridge south of Australia 
is now much better determined. 

Few measurements of crustal thickness have been made in Australia. With the 
co-operation of the Ministry of Supply, observations have been made on two 
series of atomic explosions at Maralinga. The first set of observations (Bolt, 
Doyle & Sutton 1958) over a 13° straight line traverse, yielded a depth of 35 km 
to the Mohorovitié discontinuity and no evidence of layering above this level. 
Similar results (as yet unpublished) have been obtained on a second traverse. A 
depth of the order of 36km has also been found from observations of a number 
of large quarry blasts at the Adaminaby dam site in the Snowy Mountains region 
(Doyle, Everingham & Hogan 1959). 

Some seismic work has recently been done at sea off the southwestern and 
southern coasts of Australia by the Lamont Geological Observatory’s Research 
ship Vema in association with one of the Royal Australian Navy’s hydrographic 
vessels HMAS Diamantina. Four reverse traverses were made on the continental 
margin off-shore between Fremantle and Cape Naturaliste, W.A.; four reverse 
traverses at longitude 116° E from the coast to latitude 37°S and three reverse 
traverses on the continental margin south-west of Adelaide, S.A. Refractions 
from the Mohoroviéié discontinuity appear to have been recorded on three or four 
of the traverses but the results have not yet been computed. 

Work on theoretical seismology has been almost entirely confined to Professor 
Bullen and his group in the University of Sydney whose work on travel-time 
curves and the interior of the Earth are too well known to need mention (cf. Bullen 
1953, 1956a, b; Bolt 1959). Mention should, however, be made of the work of 
Burke-Gaffney & Bullen (1958) on observations of the 1954 Bikini explosion. 


Rock magnetism 


Several Australian universities are active in this subject: the Universities of 
Sydney (Narain & Bhaskara 1957), Adelaide, the University of N.S.W. (Curnow & 
Parry 1955) and the Australian National University. The Australian National 
University operates two astatic magnetometers, one suitable for igneous and one 
for sedimentary rocks, the latter with a period of 20s and a sensitivity of 6 x 107 
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mm/gauss. A spinner magnetometer with a platinum-rhodium specimen holder 
operating at 78cps in a furnace at temperatures of up to 950°C has also been 
built; this is capable of measuring moments of 5 x 10 e.m.u/cm® (Stacey 1959). 
There is also a susceptibility balance for measuring Curie points, and apparatus for 
demagnetizing rock specimens by alternating magnetic fields and by heating. 


Table 2 


Availability in Australia of rocks suitable for palaeomagnetic observations, 

that is, igneous rocks and red sediments. 1: readily available; 2: available 

without difficulty; 3: difficult; E: rocks extensively sampled; E?: date 

uncertain; P: preliminary sampling only has been made. Measurements 
are in progress on those rocks listed as having been sampled. 


Sampling Availability Reference 


Quaternary : Irving & Green (1957). 
Green & Irving (1958). 

Tertiary 

Cretaceous 

Jurassic OF Irving (1956). 

Triassic : Irving & Green (1958). 

Permian > 

Carboniferous D I Irving & Green (1958). 
Irving (1957 a, b). 

Devonian » Irving & Green (1958). 

Silurian Irving & Green (1958). 

Ordovician 

Cambrian : Irving & Green (1958). 

Proterozoic E 


The major activity to date has been a stratigraphic survey to determine the 
direction of magnetization of magnetically stable, well-dated rock formations and 
to deduce the position of the axis of the Earth’s dipole field relative to Australia 
at the corresponding time and to compare this with equivalent results from other 
continents. ‘The availability in Australia of suitable well-dated rock formations 
for such a survey and its progress are set out in Table 2. The results may be 
expressed either in the conventional form of a plot of the movement of the mag- 
netic pole relative to the present distribution of continents or as a plot of the 
variation with time of the latitude of the continent and of the angle which a line on 
it drawn in the present NS direction made with the NS direction at any earlier time. 
The former plot has been given by Irving & Green (1958a) and the latter is shown 
in Figure 3 below. 

In this survey repeated reversals have been found to be common through a 
wide span of geological time notably in the Pre-Cambrian and the Cainozoic. 
However, at other periods the polarity is constant, for instance the whole of the 
Permian appears to be reversed. 

Many rocks have been found in which a secondary magnetization along the 
present Earth’s field masks their primary magnetization, the primary magnetiza- 
tion being that acquired at or about the time of formation. Detailed studies of the 
effect of alternating fields have shown that in most volcanic rocks the unwanted 
secondary magnetization can be removed by treatment in alternating (peak) fields 





Geophysics in Australia 457 


of 100-300 oersted and this technique is of great assistance in the stratigraphic 
survey mentioned above. 

The fact that the Mesozoic dolerites of Tasmania are magnetized almost 
vertically makes bore core from them particularly useful for a study of the varia- 
tion of the magnetization within a rock mass. It has been used in this way by 
Jaeger & Green (1956), Almond, Clegg & Jaeger (1956) and Jaeger & Joplin (1955). 
The last-mentioned paper gave evidence for the existence of a reversal within a 
single intrusion, but this has not been observed in a hole recently drilled 500 ft 
from the earlier one. Green & MacDougal (1958) have used magnetic measure- 
ments to study surface rock movements which are of interest to engineers; for 
example, they were able to state that the rocks in the first 35 feet of a dolerite bore 
core had undergone irregular tilts averaging 4° but that below this level no move- 
ment had taken place. 
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Fic. 3.—The variation in latitude of a point in S.E. Australia (Wagga, 
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ments on Australian rocks. The arrows show the former directions of a 

line in the continent drawn in the present NS direction; the widths of 

their heads correspond to 95 per cent confidence limits (this representa- 

tion becomes unsatisfactory for high latitudes, so arrows for such points 
are not shown). 


In addition to the work on Australian rock formations, work has commenced 
on intrusive and sedimentary rocks from Antarctica (Bull & Irving 1960a, b) work 
on which is continuing, and on the volcanic rocks of Aden and the Pacific Islands. 
The results obtained from Antarctica together with those obtained by other 
workers give a preliminary polar curve relative to Antarctica which is different 
from that for Australia and other continents. 

Studies of thermoremanent magnetization of igneous rock samples induced 
in the laboratory under stresses of up to 1000kg/cm? have been made (Stott & 
Stacey 1960), and it has been found that the direction of magnetization coincides 
with the field direction and is independent of the stress. A torque-meter method 
of measuring the magnetic anisotropy of natural rocks has been developed (Stacey 
1960) which allows alignment of magnetic grains of less than 1 part in 10° to be 
measured. The alignment of crystallographic axes of cubic crystals as well as of 
grain shapes can be determined by Fourier analysis of torque curves. 
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Heat flow 


Few measurements of this quantity have been made in Australia until quite 
recently. Newstead & Beck (1953) logged all available boreholes in Tasmania and 
obtained relatively high values for the geothermal flux ranging from 2 to 2°5 u 
cal/em*s. With the collaboration of the Snowy Mountains Authority a very 
extensive series of measurements of temperatures in boreholes and tunnels in the 
Snowy Mountains region is being made. Measurements in seventeen boreholes in 
the region give values scattered about a mean of 2-2.cal/em?s (Beck 1956) and a 
figure of the same order has been obtained from a preliminary reduction of tempera- 
tures in one tunnel; in both cases these results are, as yet, uncorrected for uplift. 

Recently, Mr L. E. Howard has made an extensive series of measurements in 
mines and drill holes at a number of sites throughout Australia. The results which 
are now being prepared for publication show “normal” values of the heat flux 
with a mean of around 1-2 cal/cm?s at seven sites in the Western Australian portion 
of the continental shield while thirteen sites in other parts of Australia give rela- 
tively high values with a mean of about 2-1 ucal/cm?s. 

A great many measurements of rock conductivities have been made, some 
with a divided bar apparatus or a modification of it (Beck 1957) and some by 
transient methods (Jaeger 1959a). A study has also been made of the possibility 
of measuring rock conductivities im situ, in particular, in water-filled boreholes 
(Beck, Jaeger & Newstead 1956). 


Isotope abundance and age determination investigations 


Studies of the isotope abundances of elements and the determination of the 
radioactive ages of rocks and minerals of geologic interest have been carried out in 
the Physics Department of the University of Western Australia, using a variety 
of mass spectrometric techniques. 

This work commenced with an examination of the relative abundances of the 
stable carbon isotopes C!2 and C}% in nature (Jeffery & others 1955). As a result of 
this original work the tentative suggestion was made that there might be a measur- 
able age effect in the C!2/C!8 ratio in carbon of different ages due to systematic 
alterations in this ratio in the oceanic reservoirs. Further measurements (Compston 
1960) have substantiated this theory, but very extensive sampling would be 
necessary to provide completely convincing evidence of its validity. 

Several reports of variation in the abundances of the potassium isotopes 
initiated a study of this effect which led to the conclusion that many of the dif- 
ferences reported may have been due to sample contamination. No significant 
variations in the abundances of these isotopes were found in any of the wide 
variety of samples investigated, provided only pure potassium salts were used for 
mass spectrometric measurements (Kendall 1960). 

Potassium from minerals which gave anomalous potassium-—argon ages was 
examined to determine whether an abundance difference in the radioactive parent 
isotope was causing an age peculiarity to be observed. All the minerals examined 
contained only normal potassium. 

An exploratory survey of Australian mineral ages using the uranium-lead 
method provided the basis for further investigations in particular areas of geologic 
interest (Greenhalgh & Jeffery 1959). More recently potassium—argon and rubi- 
dium-strontium ages have been measured to substantiate and extend the uranium- 
lead measurements (Wilson & others 1960). 
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A method of determining the metamorphic history of granite rocks has been 
developed leading to the possibility of measuring both the age of a rock’s emplace- 
ment and the time at which a subsequent metamorphism took place (Compston & 
Jeffery 1959, 1960). This work is being extended to other rock types and to the 
investigation of specific geologic situations. 

At the Australian National University, pure samples of zircon have been 
extracted from acid rocks obtained from the Northern Territory, and from the 
Mount Isa and Tumut areas. Preliminary estimates of age have been obtained 
by a combination of “‘lead-«’’ and chemical methods. A gas-source mass-spectro- 
meter has been applied to the study of Australian ore-leads. A convenient method 
has been developed for the production of tetramethyl lead, which works equally 
well for galenas or oxidized ores, and which operates with a minimum of attention. 
Precision of about 0-3 per cent on the 204 isotope is currently being attained. 
Principal areas studied to date are the Northern Territory and the Mount Isa 
region. It is hoped to extend the programme greatly to cover all the usual methods. 
Work on age-determination using the potassium—argon method is being begun at 
the University of Sydney. 


Deformation and phase equilibrium studies at high pressures 


Research on deformation of rocks and minerals at pressures of up to 10000 
bar is in progress at the Australian National University (Paterson 1958, 1959; 
Jaeger 1959b). A “‘squeezer”’ type apparatus is in use by A. E. Ringwood for study- 
ing phase equilibria among magnesium, nickel, cobalt and manganese silicates 
and germanates at temperatures up to 1 500°C and pressures up to 80000 atmo- 


spheres. The object of these studies is to gain further data on the olivine—spinel 
transition and its role in the mantle. The nickel analogue of olivine, NigSiO,, has 
been found to invert to a spinel structure around 30000 atmospheres at 650°C. 
Methods of calculating pressures required to cause phase transitions are also being 
developed. 


Australian National University, Bureau of Mineral Resources, 


Canberra: Melbourne. 
1960 June. 
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1960 July 26 to August 6 


Thanks to the co-operation of a number of their colleagues, the Editors present 
below reports of some of the more interesting sessions of the Assembly. With seven 
Associations holding simultaneous sessions, it is not possible to report everything, 
but we hope we have chosen some of the highlights. 


1. Symposium on the Geodetic Uses of Artificial Earth Satellites 


The first session of this symposium on July 29, under the chairmanship of 
Professor Roelofs, began with the presentation of a report by Dr W. Markowitz 
on the geometrical uses of artificial satellites. If observers at groups of stations 
determine, say, the topocentric directions of a satellite from these stations, either 
simultaneously or at known times, they can derive the relative positions of the 
stations geometrically, although if the observations are not simultaneous the orbit 
of the satellite has to be known. Other methods depend on measuring the distance 
of the satellite or its velocity relative to the observatories. Dr Markowitz suggested 
that for optical observation the light on a satellite should be intense enough to be 
photographed with a 10- to 15-in. telescope, while a flash rate of one per minute 
had been proposed. Experiments with a Xenon lamp had given flashes 1 ms 
long with energies of 32000 lumen. One would expect to determine topocentric 
directions with an accuracy of 1” to 2”. Although, as Dr Markowitz pointed out, 
no systematic geometric observations had been made, since so far no satellite in- 
tended for this purpose has been launched, Dr Whipple mentioned in the discussion 
that photographs of satellites already launched, taken with the Baker—-Nunn 
cameras, were being analysed at the Smithsonian Astrophysical Observatory to 
derive the relative station positions and later in the discussion Dr G. Veis gave 
some details of this analysis and said that the standard deviation of the computed 
positions referred to the International Ellipsoid was 50m. 

An important point brought out in the discussion was the effect of atmospheric 
refraction and scintillation. Dr R. d’E. Atkinson remarked that much greater 
accuracy than 1” was possible with an astrographic camera but that with a single 
flash at 1-min intervals the error due to scintillation would be serious and a group 
of flashes should be given to average out this error. There should be several per 
second. Dr Whipple added that the standard deviation of 2” for a single observa- 
tion was essentially due to scintillation and agreed that a group of flashes was 
wanted. He also remarked that because of the limitation imposed by scintillation 
there was no gain in using a very long focus camera. Another important considera- 
tion for optical observation was the number lost because of clouds. Dr Whipple 
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said that with the present satellites, for which observations were limited to twilight, 
on the average 12 Baker-Nunn observations were obtained per month on each of 
8 satellites while Woomera gets almost 24 observations per satellite per month. 
He estimated that with a flashing light twice or thrice as many observations could 
be made. Dr A. M. Stone said that with radar or Doppler methods there are no 
such limitations and already enormous quantities of data are being obtained from 
Transit, the navigation satellite using Doppler velocity measurement, with which 
geodetic measurements can be made to 4 or 5m. 

Dr A. H. Cook presented the second report, on the determination of the Earth’s 
gravitational potential from observations of artificial satellites, in which he sum- 
marized the principles on which the external potential may be derived from changes 
in the orbits of satellites. He gave an outline of the methods and results of theories 
of the perturbations of orbits by the Earth’s gravitational field and mentioned 
theoretical problems still outstanding, and summarized the principal numerical 
results obtained. He pointed out that the estimates of the coefficients of zonal 
harmonics in the potential were about twenty times as accurate as those derived 
from harmonic analysis of surface gravity measurements but emphasized that values 
of certain of the higher ones could not be accepted with complete confidence until 
certain discrepancies between theories had been resolved. Among a number of 
important points made by Dr W. M. Kaula in the discussion was that the value of 
the third harmonic from Transit observations, reduced on a purely numerical basis, 
confirmed the value derived from Vanguard I data and so checked the theory on 
which the latter depends. He also pointed out that the real limitation to more 
accurate analyses of orbital changes was the irregularity in atmospheric drag. 
The geophysical significance of some of these results was also mentioned in the dis- 
cussion, in particular the demonstration that the state of stress in the interior of the 
Earth was not hydrostatic. 

In the remaining part of the first session of the Symposium, two papers were 
discussed—one by I. D. Zhongolovitch, presented by Professor A. Mikhailov, 
giving particulars of his theoretical work and some new numerical results, and one 
by Dr E. Buchar who has applied Zhongolovitch’s theory to the data from five 
satellites. 

The second session of the Symposium on August 1 was devoted to papers on 
the use of satellites for making geometrical connections. Dr M. Henriksen discussed 
some studies of results from Minitrack Mk II tracking Vanguard I in which the 
longitudes of four stations had been determined with an uncertainty of between 
100 and 500m with respect to N. American datum using observations lasting two 
to three months. The aerial patterns of the Minitracks had been calibrated by 
observations of radio stars. Higher accuracy will be obtained when the radio 
frequency is 216 Mc/s instead of 108 Mc/s as at present, but errors due to un- 
certainties in the determination of the orbit will still amount to about + 200m. 
In the discussion of this paper Dr A. M. Stone gave some particulars of the 
Doppler velocity measurement system used in Transit. The accuracy of results 
so far obtained is about + 40m for a single passage, not as good as the ultimate 
possible but considerably better than Minitrack. The reason for these good results 
is that the quartz crystal oscillator carried in the satellite has the remarkably good 
stability of about 1 in 109 over short periods. ‘Two frequencies, 216 and 324 Mc/s, 
are used so that ionospheric refraction errors can be eliminated. 

The optimum use of a geodetic satellite was considered by O. W. Williams 
who discussed the distributions of stations, the limitation on optical observations 


2c* 





464 XII General assembly of the International Union of Geodesy and Geophysics 


and the advantages of combining them with electronic velocity and range measure- 
ments and the desirable properties of orbits. 

Dr Hertz presented a paper in which he reviewed some applications of celestial 
mechanics to satellite problems, mentioning the importance of defining precisely 
the relation between observations and the parameters that enter theoretical work. 

Following discussions at this symposium and otherwise, the International 
Association of Geodesy has set up a Permanent Commission on the geodetic uses 
of artificial satellites, under the presidency of Dr A. H. Cook, with the objects of 
advancing theoretical studies and the organization of geodetic observations. 


2. Symposium on Island Arcs 


This symposium was held under the chairmanship of Professor Vening Meinesz, 
who received an ovation on opening the meeting. 

Professor V. V. Beloussov presented a paper by himself and E. M. Rounditch 
on “The position of Island Arcs in the development of the Earth’s structure” in 
which they discuss the form of island arcs in the light of the hypothesis that the 
Earth began as a cold agglomeration and has warmed up through radioactive heat- 
ing, and that the oceans are comparatively recent features that have formed through 
the basification of the continental crust and its replacement by oceanic type of crust. 
Island arcs could be classified into two groups. Those in the first are similar to 
folded mountain ranges and were formed from geosynclines, but in the process of 
basification of the crust, the interior massifs of the geosynclines have subsided, 
leaving arc-like relics of continental crust, such as Japan, Indonesia, or the Antilles. 

The second group comprises island arcs with no geosynclinal antecedents and 
were formed as the direct result of crustal tension. They are younger than arcs of 
the first type. The curvature is due to the distortion of the stress trajectories in 
the crust due to the initial fracture of mantle material. The fracture facilitated 
differentiation and the rise of light material to the surface. Examples of such arcs 
are the Aleutians, the Kurile Island and the Marianas. Although the differentia- 
tion of mantle material has on occasion produced surface results similar to geo- 
syncline development, island arcs cannot be identified in the geosynclines. Profes- 
sor Beloussov pointed out the asymmetry of the borders of the Pacific Ocean, the 
Eastern being a single zone with a unified development while on the West, the crust 
has a mosaic structure with more complex development. 

In reply to Professor Vening Meinesz, Professor Beloussov thought that convec- 
tion currents occurred in the mantle but were not so much the result of thermal 
forces as a manifestation of chemical differentiation. 

Dr J. Hennion presented a paper by himself and J. and M. Ewing on “‘Seismic 
refraction studies of island arcs and oceanic trenches’’, dealing especially with the 
West Indies arc, the Puerto Rico trench and the Cayman trough. The island arcs are 
characterized by great thickness of crustal rocks with a longitudinal seismic velocity 
increasing with depth from 4 to 8km/s. The Mohoroviéié discontinuity is depressed 
below both island arcs and trenches, to about 40km below the former and 20km 
below the latter. The thickness of sediments varies from some kilometres in the 
Puerto Rico trench to almost nothing in the South Sandwich and Tonga trenches. 
Taken with these seismic results, gravity anomalies given in the paper ‘‘Gravi- 
metric studies of ocean trenches” by M. Talwani, J. L. Worzel and M. Ewing, show 
that the great variations of anomaly are predominantly due to the variation of the 
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depth of the sea and of the Mohorovitié discontinuity, variations in sediment and 
crustal densities having only minor effects. As compared with isostatically com- 
pensated structures, there is an excess of mass under Puerto Rico and a deficit 
under the trench. Similar comparisons of seismic results and gravity anomalies in 
the Tonga—Kermadec region give rise to difficulties in fitting suitable structures 
although one possibility is that there is a low velocity layer above the Moho so that 
the apparent depth of the Moho is too small. 

Professor Glangeaud, commenting on the previous papers, pointed out that 
both types of arc appear to occur in the Mediterranean. 

Russian results in the Kurile region were presented in a paper ‘““The structure 
of the Earth’s crust in the region of the Kurile Island arc’, by Messrs. I. P. 
Kosminskyn, K. M. Krakshinn, G. G. Mikliofa and J. V. Tulina. They show that 
crustal structure can be divided into three types according to clear distinctions 
between seismic travel-time curves and velocity depth curves in these types. On 
this basis, the Kurile Island arc can be divided into three parts in each of which 
the crustal structure is similar to the adjacent oceanic or continental region—thus, 
in the North, Kamchatka is of continental type. There are corresponding dif- 
ferences in gravity anomalies and in the depth of the Moho. Earthquakes occur 
mostly in the continental region while volcanic activity is greatest in the inter- 
mediate area. 

Further work on the South Kurile Islands was presented by S. A. Fedotor in a 
paper “Preliminary results of a detailed seismic investigation of the Southern 
Kurile Islands”, by himself, V. N. Averianorn, A. M. Bagdasaror and I. P. Kurzin. 
The seismicity of this region was studied with a special network of five high 
sensitivity seismic stations having a maximum response at 5 c/s. Co-ordinates of 
epicentres were determined to within 1okm. It has been found that for epicentres 
at depths of from 30 to 45km the attenuation constant is 0-0047km~-! while for 
the range 30-85 km it is 0-o20km~!, or about 0-03 for the period of maximum 
response (5c/s). These values indicate that the material shows more non-elastic 
behaviour with increasing depth. The epicentres are scattered about a plane 
sloping down at about 45° to the west, with the epicentres confined to the region 
below the strip between the Kurile arc and the deep trench. Only shallow focus 
earthquakes occur in the volcanic region of the Kurile arc; apart from this, all 
epicentres lie below the Moho. 

An account of morphological studies was given by R. W. Fairbridge in his 
paper, ““The Melanesian border plateau”’, a zone of crustal shearing in the S.W. 
Pacific. He called attention to the fact, indicated at the Toronto assembly, that the 
majority of fault displacements around the Pacific turn the basin in an anti- 
clockwise direction. 


3- Mechanical Constitution of the Earth 


3.1. Structure of the Earth’s crust 

Markus Bath has made a useful summary of all available seismic results relating 
to the Conrad discontinuity in continental rocks. He concludes that there is a 
general increase in velocity downwards in the crust and that a Conrad layer does 
seem to exist. The divergences in the results are almost exclusively due to the 
interpretation of the records and the travel-time plots, and the observations do not 
justify the conclusion that different continents have different structures. In a 
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discussion of these remarks it became clear that part of the trouble with conven- 
tional continental seismic observations is that not enough use is made of any waves 
other than the first arrival. As Hales pointed out, the need for the use of closely 
spaced arrays in order to correlate later arrivals has been recognized in reflection 
work, but except for the work of Gamburtsev (whose book Correlation Methods in 
Refraction Work has recently been translated in full into English and should be 
read by all those who do refraction work) this has not been so with most refraction 
experiments, and in very few cases has the amplitude been used as a guide to the 
identity of waves. Willmore added another word of advice to those embarking on 
long-distance refraction work. He underlined the fact that experiments should be 
so laid out that the relative value of the intercept due to the travel time under 
each explosion and geophone is discovered. It is important to follow the pre-war 
teachings of Bullard and to ensure that there are enough geophones with close 
grouping around any location, so that a good enough velocity measurement can be 
obtai.2d to show whether variations from shot to shot are due to refracted waves 
going deeper or to sideways variations in depth of refractor. 

Steinhardt, Meyer and Woollard contended that much of the ambiguity in the 
results of refraction surveys can be removed by the proper use of the times of 
arrival and amplitudes of reflected waves, supplemented by gravity anomalies. They 
illustrated their paper by discussing the results of surveys in Mexico, Wisconsin 
and Arkansas. 

Frank Press showed how surface wave phase velocity observations may with 
advantage be combined with refraction and gravity work. Checks on the mean 
shear wave velocity and on rock density in the crust may be obtained at places where 
good refraction results are available and the surface waves may be used to extend 
the “rustal thickness measurements cheaply and rapidly over a large area. Press 
and Knopoff are at present engaged in setting up surface wave observation stations 
in Europe. The results reported for the oceans suggest that the Pacific, apart from 
the Easter Island rise, is slightly different in crustal structure from the Atlantic. 

A great deal of work on thickness of the Earth’s crust has taken place in the 
U.S.S.R. during the past years and it was of great interest to hear from G. K. 
Tvaltvadze that the original method used by Mohoroviéié—of determining crustal 
thickness from observations of near earthquakes—gives results which are in fair 
agreement with the results obtained using artificial explosions and with thicknesses 
estimated from gravity anomalies. The virtue of near earthquake observations is 
that small earthquakes occur mainly in mountainous regions which are just those 
areas which offer little opportunity for the other methods of investigation. 

A different approach to the use of earthquake waves was discussed by A. A. 
Treskov, who has made extensive measurements of crustal thickness by observing 
waves from deep-focus earthquakes reflected from the base of the crust. 

Professor G. P. Woollard described an estimate of crustal structure in Antarc- 
tica based on the seismic measurements of ice-cap thickness combined with his 
hypothesis (backed by many comparisons of gravity and seismic structures all over 
the world) that most of the structures in the world are in isostatic equilibrium and 
that the density distribution with depth is approximately constant. The I.G.Y. 
gravity measurements treated on this basis suggest that although much of Western 
Antarctica was an island archipelago prior to being covered by the present ice-cap, 
the area as a whole is characterized by a continental type of crust varying in thick- 
ness from about 28km to over 40km.* Assuming that removal of the ice-load 

* But see Evison & others (1960) Geophys. J., 3, 289. 
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would entail a readjustment to maintain isostatic balance, it is calculated that a rise 
in the land mass of about o-gkm would take place. 

Of the results described, the surveys in the Russian platform and Central 
Asia presented by Y. N. Godin appeared to have assimilated all the advice from 
earlier papers and also from Gamburtsev. Geophone spacing was kept down to 
100m, while the distance between shot points was 40-80km. Altogether 3000km 
of crustal structure investigation have been studied in the last three years. The 
crustal studies in the Alps (sub-committee for the European Seismological Com- 
mission on Alpine Explosions) and in the Rocky Mountains (Meyer, Steinhardt, 
Woollard & Bonini) are not yet complete, but it is clear from the results obtained 
to date that there is no indication of an appreciable increase in the depth of the 
Mohorovitié discontinuity under the Rocky Mountains, whereas there is a distinct 
dip towards the centre of the Alps. 

In Iceland, where Bath made seismic refraction measurements, the crustal 
structure is intermediate between continental and oceanic. 

The dispersion of surface waves has also been used to find crustal structure by 
Oliver, Sutton and Ewing, who illustrated the various types of seismogram that 
may be obtained by presenting records from some of the stations co-operating in 
the I.G.Y. long-period seismic programme. By analysing the dispersion curves 
for the various modes, including Rayleigh waves and Love waves, with periods 
greater than 14s, corresponding to fundamental modes, and other surface waves 
and P waves probably guided by the crustal layers, they were able to obtain 
a good idea of the structure along the path of the waves. Ocal and Sikharulidze 
described the determination of the structure of various parts of Asia using this 
technique. 

Long period (30~74s) Love surface wave results were also discussed by I. I. 
Popov, who compared waves along the shortest path with those travelling the 
opposite way round the Earth. The results agreed with crustal thicknesses of 
35 km under the continents, 15 km in intermediate areas and 6km under the oceans. 

A complicated crustal structure was reported by A. G. Averianov for the 
northwest corner of the Pacific. ‘There seems to be an increasing volume of seismic 
evidence for the divisions of ocean crustal structure to include an intermediate 
type, as proposed by Gaskell. Both this Russian work and that of a combined 
4-ship team in the Caribbean were in the opinion of many patient listeners, repeated 
too many times by different speakers. A good summary of several papers in both 
instances would have made an interesting and understandable exposition and would 
have allowed time for questioning and discussion. 

The results of the Caribbean experiment have been described under Island 
arcs. Similar seismic measurements by E. Bunce and D. Fahlquist of Woods Hole 
Oceanographic Institution were the first to find the depth of the Moho in the 
Mediterranean, where Gaskell and Swallow and Ewing & others had previously 
failed to get below a 4-5km/s layer. ‘The structure found in the Tyrennian Sea 
shows a good velocity of 7-7km/s beneath 6-6km/s and sediment plus 4—5 km/s 
material. The low mantle velocity could well be the result of the past history of old 
sedimentation followed by continental drift proposed by Glangeaud. 

Bott has made a study of gravity anomalies over granite masses and finds a 
consistent pattern of 50-80 mgal negative anomalies. His computer method of inter- 
preting the observations suggests a structure above the Conrad discontinuity of 
10km of granite beneath a further 1okm metamorphosed sedimentary rocks. 
The granite has a lower density than the metamorphics, and granite masses are 
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manifestations of the lighter granite rock forcing itself up to the surface through the 
overlying sediments. 

In reporting correlation of climatic zones in cores from the deep ocean floor, 
D. Ericsson stressed the need for much longer cores if the geological history of the 
oceans is to be unravelled. Good news in this respect was given by Dr A. Maxwell 
who announced that the AMSOC Committee, with financial support from the 
National Research Council, planned in January 1961 to drill about 500 feet into 
the sediment at a location near Guadeloupe Island, where the water depth is 
12000 feet. ‘The oil company drilling barge Cuss J will be used for this operation, 
the purpose of which apart from obtaining information about deep sediments, is 
to discover the difficulties associated with drilling from an unanchored ship with 
several miles of drill-pipe hanging free in the water. 

Although the Mohole project is the implementation of a resolution made at the 
XIth Congress at Toronto in 1957, it also forms part of a new resolution passed at 
Helsinki. This is a project sponsored by Professor Beloussov: ‘The Upper Mantle 
and its influence on the development of the Earth’s Crust.”” The Russian interest 
in putting forward this broad target for seismological and oceanographical research 
may, it is hoped, indicate that funds are to be made available for drilling into the 
Earth’s mantle. In the opinion of many of those associated with oil-well drilling 
the successful penetration of the Earth’s crust is merely a matter of money so that 
some definite checks of the geophysical results may be achieved before the next 
1.U.G.G. Congress. 

The ocean sediments are being investigated with new and improved geophysical 
tools. Dr J. B. Hersey has obtained many traverses of topographical depressions in 
the oceans with the precision echo-sounder. This instrument measures changes of 
a few feet in the depth of the 3-mile deep ocean, and often shows sub-bottom echoes 
which show stratification of more compacted sediment layers. The abyssal plains 
have been classified roughly into large plains which face the continents and which 
show no consistent multiple sub-bottom echoes extending more than a few miles 
along the ship’s track, and small, totally enclosed basins ranging from a few miles 
up to 100 miles in size, which show continuous layering of sediment across the 
basin. The Woods Hole research vessel Chain was on view at Helsinki, complete 
with the latest versions of the “Sparker’” and ‘‘Exploder’” devices for deeper 
penetration of sea-bed sediments. This (the ““Thumper’’) is a towed diaphragm 
type of no’se-maker, whose power output approaches that used in conventional 
seismic reflection work, but which can obtain continuous profiles by sending out 
sound impulses every second. 

An example of the enormous quantity of facts that are being collected over 
the years by the indefatigable Professor Maurice Ewing and his school at Lamont 
was given by Dr J. Worzel. By combining continuous gravity profiles with seismic 
observations, it is shown that, on the eastern seaboard of the U.S.A., the true 
continental edge is the thousand fathom curve. Continental crustal thicknesses 
thin to the normal deep ocean figure in a distance of about two hundred kilometres. 
Magnetic measurements support this picture, and show, for example, the existence 
of a buried sea-mount with 4o0o0oft of relief under the continental shelf off New 
Jersey. 

Professor Ewing, as befits the pioneer of deep-sea seismic work, has produced 
yet another advance by constructing a geophone which operates on the sea-bed by 
signalling results to a surface ship by means of a modulated 12ke acoustic wave 
train. There is a good harvest of results to be gathered regarding deep-sea 
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sediments and Layer 2 using this new tool, and the geophone may possibly be 
adapted to record earthquake waves in those parts of the ocean where no convenient 
islands exist on which to set up a seismological station. 

Seismic and gravity results over the Mid-Atlantic Ridge still leave some doubt 
as to the geological structure of this feature. The basement iayer has the velocity 
7:3km/s and the gravity anomalies possibly require a low density root to satisfy 
them. 


3.2. The internal constitution of the Earth 


The late Professor Gutenberg would have been delighted to see paper after 
paper presented to the Seismological Association confirming his low velocity 
layer. It was at the 1.U.G.G. Congress in Rome that Professor Gutenberg ex- 
plained his conviction that amplitudes of P waves required a reversal of the normal 
increase of velocity with depth in the mantle. The new evidence provided by the 
dispersion of surface waves, which has been forthcoming now that several long 
period seismographs are in operation, fully supports the concept of a zone of low 
shear velocity between 100 and 200km beneath the continents. Under the oceans 
the region of low shear velocity rises to about 50km depth, and the low velocity 
zone is thicker, and the velocity lower, than under continents. Careful examination 
of travel-time curves by Miss I. Lehman confirms this picture, and the weight of 
the evidence suggests that the boundaries of the low velocity layer are sharp 
discontinuities in velocity, and probably in density. 


3.3. Periods of free oscillations of the Earth 

In 1953 Benioff reported the possibility that he had observed a free oscillation 
of the Earth on his strain seismograph recording of the Kamchatka earthquake of 
1952. This announcement stimulated many theoretical calculations—by Pekeris, 
Stoneley, Jobert, Takeuchi, and Gilbert and Macdonald among others. 

In May 1960 Press read a paper to the American Geophysical Union reporting 
an analysis of the Press-Ewing vertical seismogram (79 = 15s Ty = gos) from the 
Kamchatka earthquake. He had associated the following spectral peaks with 
higher mode spheroidal oscillations: 9-65, 8-34, 7:46, 6°49, 5:26min, but he con- 
sidered that verification from another earthquake was needed. Bogert supplied 
this verification at Helsinki when he announced the result of his analysis of the 
recordings of the great Chile earthquake of 1960 May 22 by a similar instrument in 
New Jersey. The periods from the Kamchatka and Chile records were almost 
identical. 

The most exciting news on this subject came from Benioff (by telegram) and 
Slichter, who had recorded the fundamental spheroidal mode and a series of low 
modes from the Chile earthquake. Benioff and Smith’s analysis of the Benioff 
strain seismogram from Isabella (California) had revealed surface waves lasting for 
two days after the time of the shock. The power spectrum of these waves showed a 
long series of peaks: the greatest period observed was 53-1 minutes, and this was 
followed by periods of 35-7, 25-8, 21-8, 19°8 minutes, and in regular sequence 
down to 6-01 minutes. These observations were found to be in remarkable agree- 
ment with the theoretical periods forecast by Pekeris in his analysis of the free 
vibrations of an Earth of the constitution of Bullen’s Model B. Four torsional 
modes were identified, and 15 spheroidal. 

Slichter reported that the gravity meter of high sensitivity specially constructed 
for the same , urpose at U.C.L.A. had also recorded long period oscillations due to 
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the Chile earthquake. No torsional modes had been observed, but the periods 
of spheroidal modes agreed to 0-1 min with those reported by Benioff, Press and 
Smith.* There is little doubt that the new experimental facts, for example the 
observed doubling of the spectral lines, will inspire theoretical studies to move 
another step towards a true understanding of the Earth’s internal structure. 


Table 1. 
Results of power spectra analysis, Chilean earthquake 1 


Observed periods in minutes 
Spheroidal motion Toroidal motion 
(absent on Z) 


54°7 7*10 42°3 
53°1 6-78 28:6 
35°9 6°49 218 
35°2 6°23 17°9 
25°8 6-01 15°4 
24°5 5°78 13°5 
19°8 5°59 12°3 
14°0 5°39 11-2 
(151) 5°26 (11-0) 
12°5 5‘10 10°3 
11°8 4°96 
10°6 4°83 al me 
(9:8) 4°69 There are additional peaks on 
9°6 4°59 the spectra with high coherence 
8-94 4°47 between Nana and Isabella: 
8°37 4°27 75+3 
7°88 4°18 100 +5 
7°47 4°08 
Instruments: Isabella strain NW 
Nana strain N-S 
Pasadena 30-90 UD 


Pasadena 30-90 NS 
Pasadena Exper. EW 


Pekeris pointed out that the gravimeter would not respond to torsional modes, 
so that the agreement with his theoretical results was excellent. He further outlined 
theoretical methods of obtaining the periods of all Rayleigh modes, from the 
fundamental down to those recorded by ordinary seismographs, when a particular 


model of the Earth supplies the necessary hypothesis on variation of elastic para- 
meters with depth. 


4. Age Determination in Rocks and Meteorites 
A series of meetings on aspects of geochemistry was arranged at Helsinki in 
conjunction with the International Associations of Seismology and Vulcanology. 
Besides the presentation of a number of scientific papers, the meetings aimed at 
convening a representative group of geochemists to discuss the international 
organization of this science. This report summarizes some interesting aspects of 
the meetings on age determination and meteorites. 


* Hugo Benioff, F. Press & S. Smith, ‘‘Experimental verification of the free oscillation of the 
Earth’’, J. Geophys. Res. 1961, February. 
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4.1. General 


A number of laboratories are now equipped to obtain several figures for the 
age of a rock by using two or more radioactive decay schemes on various minerals 
separated from the rock. Typically a rock may give Rubidium-Strontium ages 
on feldspars and micas, Potassium-—Argon ages on micas, and U-Th-Pb dates on 
zircons. It is rare to obtain the same figure from all the independent determina- 
tions. L. T. Aldrich (Washington, D.C.) gave some examples which showed how 
varied are the patterns of ages which may occur in different geological settings. He 
suggested the following order of resistance to change, beginning with the method 
which appears to give the most stable value of the age. 


. Feldspar (Rb-Sr) 6. Feldspar (K—A) 

. Muscovite (Rb-Sr) 7. Zircon (U235—pp207) 

. Muscovite (K—A) 8. Zircon (U238—Ph2%6) 

. Biotite (Rb-Sr) 9. Monazite (Th?32—Pb298) 
. Biotite (K-A) 10. Zircon (Th232—Pb2%8) 


It is not uncommon to find the time of two widely separated metamorphic episodes 
recorded fairly accurately in the same rock. However, in such a case supporting 
evidence is necessary in order to arrive at convincing conclusions, since ages inter- 
mediate between those of the two metamorphic episodes may otherwise be deduced. 


4.2. Zircons 


Measurements on zircons yield a lead ratio age (Pb2°7/Pb2°6) in addition to those 
shown above. Although this is not independent of the two U-Pb ages it is the least 
subject to change, as A. J. Burger (Pretoria) and Aldrich showed. ‘The Thorium— 
Lead age is almost invariably the worst. Loss of radiogenic lead during the life of 
the mineral is the probable cause of the discrepancies but the mechanism of loss 
is not understood. The pattern of ages which would be produced by continuous 
diffusion or leaching is not very different from that which would be produced by 
loss during a metamorphic episode. There is evidence in support of both 
mechanisms. 

S. Deutsch and C. T. Silver (Cal. Tech.) described experiments in which the 
outer layers of zircon crystals were “peeled off” by partial meltings. The ages 
from the stripped layers differed from each other and from the age of the core, but 
the pattern which they formed fitted well with the pattern obtained from the 
different size fractions of the original zircon crystals and from the associated small 
quantity of Uranothorite. This pattern could be explained by metamorphism 
80 m.y. ago of material originally crystallized 1650 m.y. ago, which is in agree- 
ment with field evidence. 


4.3. Retention of Argon 

Fechtig and Zahringer (Heidelberg) described their experiments on the diffu- 
sion of argon-37 out of calcium containing minerals in which it was produced by 
neutron irradiation. The sensitivity of the method is so high that measurements of 
diffusion at room temperature can be made. The relationship between diffusion 
coefficient and temperature was studied for several materials, but only for chondrites 
did log D vary linearly with (O/RT). For these materials losses at room tempera- 
ture are estimated to be very small. For rock-forming minerals such as augite, 
anorthite, and calcium mica, the form of the relationship is not fitted by any simple 
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mathematical formula, and it is not yet possible to explain it qualitatively. How- 
ever, these diffusion measurements agreed roughly with the retentivity found in age 
determination. 

Important results on sanidines were presented by Baadsgaard and his col- 
leagues. Potassium argon ages were compared for several pairs of biotites and 
sanidines obtained from bentonites (fine-grained clay-like rocks believed to be 
derived from volcanic detritus), and for all pairs the two ages agreed within the 
limits of error (5 per cent). Furthermore, experiments on one of these sanidines at 
an elevated temperature gave an activation energy for the diffusion process com- 
parable with that found in micas. On the other hand, evidence from a trachytic 
lava flow showed much greater diffusion losses in a similar experiment. 


4.4. Pre-Cambrian geochronology 


Absolute dating methods are an indispensable tool for the study of the pre- 
Cambrian. Most work has been done in North America, and a rudimentary 
pattern has emerged from this. This was summarized by R. D. Russell and B. J. 
Giletti. The oldest terrain, 2 500 to 2700m.y., has been found in three separate 
regions two of which are in the Canadian shield, and one in the Northern Rocky 
Mountains. Other metamorphisms occurred 1 600 to 1 800m.y. ago and 1 400m.y. 
ago, and goo to 1 100m.y. ago. The pattern of ages does not exclude the possibility 
of continental growth round the old “‘nuclei’’. 

Extensive Soviet work on the Baltic shield embodies a very large number of 
analyses, mostly by the Potassium-—Argon method. Unfortunately the authors 
(Polkanov & Gerling) of a paper on this subject were unable to attend the meetings 
in Helsinki, so that their work is available only as an abstract. Their most notable 
results are ages of 3250 to 3590m.y. for the Kola peninsula; they group these in 
the “‘Katachaean I”’ orogenic cycle. These are the greatest ages yet found. They 
list ten other pre-Cambrian orogenic cycles, but in view of their heavy dependence 
on the Potassium—Argon method, and in the absence of details of their results, it is 
not yet possible to feel complete confidence in their conclusions. 

A discussion of the geochronology of Eastern Asia was presented by Vinogradov 
and Tugarinov. This very large area has not previously been studied by absolute 
dating methods. Results of Potassium—Argon measurements on both igneous and 
sedimentary (glauconitic) rocks were presented, together with some U-Th-Pb ages. 
In addition some “‘model”’ ages were given for common lead extracted from dolo- 
mites. In summarizing their results the authors proposed the following chronology 
for the pre-Cambrian: 


Era Tectono-Magmatic Epoch 


Katangan 
Upper Proterozoic 

Grenville 1 100 
Lower Proterozoic 


Belomorian I 900 
Archaean 


Rhodesian 2700 


Questioned on their views on the time at which the Cambrian period began, 
Tugarinov considered it to be between 550 and 650 million years ago. The un- 
certainty lies partly in the difficulty of finding a satisfactory stratigraphic definition 
of the base of the Cambrian. 
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The greatest age measured outside the U.S.S.R. is that of the Old Granite 
north of Johannesburg. Hales presented details of Rubidium—Strontium measure- 
ments on this which gave an age of 3205 + 33 million years. These were “‘total 
rock” measurements; investigation of the separated minerals shows that radiogenic 
strontium has moved between them since the granite was emplaced. 

Moorbath presented work done at Oxford on the Scottish pre-Cambrian. The 
two main periods of metamorphism in the Lewisian Complex are the Scourian 
and the Laxfordian. The Scourian date has been obliterated in many rocks and 
minerals, but survives in a few. It is about 2 500 million years. The Laxfordian is 
about 1 500 million years old. The main metamorphism of the Mione and Dal- 
radian sediments is Caledonian, but there is some evidence for an earlier meta- 
morphism which may have occurred 740 million years ago. 

The question of correlation of pre-Cambrian metamorphic—orogenic episodes 
arose in discussion. Aldrich thought that the only epochs found convincingly on 
every continent are the Rhodesian (2 500-2700) and the Grenville (goo—1 100). 
However, most of the others appear in at least two widely separated regions of the 
Earth’s surface. Goldich mentioned that a valuable summary of age data has been 
assembled by Gastil (Am. Journ. Sci., May 1960) who gives histograms of the 
numbers of localities whose ages fall in each 50 million year interval of time. 
These measurements show clearly the concentration of ages in groups correspond- 
ing presumably to the growth and decay of orogenic activity in a roughly periodic 
fashion. 


4.5. New dating methods 


Two new dating methods were discussed at the meetings. One is concerned 
with the absolute dating of ocean sediments, and originates from Emiliant’s group 
at Miami. The method assumes that Thorium-230 (lonium) and ‘Thorium-231 
(Protoactinium) are precipited in identical fashion from sea water and are incor- 
porated in sediments in a ratio which depends only on the U?35/U288 ratio. It 
also assumes that the parent Uranium isotopes are not appreciably precipitated. 
The half lives of Th?8° and Th?! are 8 x 104 and 3-5 x 104 years respectively; it 
may be expected that a measure of their ratio in sedimentary material up to 
250000 years old may be used to date that material. The first results are promising. 
It is suggested that the method may be far more effective than the old ionium 
method, which depends on assumptions about the rate of sedimentation or the rate 
of precipitation of ionium which are most unlikely to be valid. 

W. Broecker (Lamont) gave preliminary results on the dating of carbonate 
materials by a new method. In carbonates uranium is taken up for more than 
thorium. Consequently after the material is formed the growth of the Th23°/U288 
ratio towards its equilibrium value should be a function of the age of the material. 
Preliminary results are promising, but not conclusive. 


4.6. Rare gases in meteorites 


In the last few years considerable work has been done on the rare gases in 
meteorites. Most of this work has been concerned with gases produced by spal- 
lation of the nucleii of meteoritic material by cosmic rays. However, there has 
recently been evidence from several workers that some rare gas isotopes must have 
been occluded in the meteorite at the time that it was formed. A paper by 
Zahringer & Geutner described some further evidence for this, including the 





474 XII General assembly of the International Union of Geodesy and Geophysics 





discovery of neon and argon with istopic ratios which could not be explained 
either by nuclear spallation or by contamination from atmospheric neon and argon. 

Important work on helium was described by Nier and his colleagues (Minne- 
sota), who measured the distribution of Helium-3 and Helium-4 contents as a 
function of position in several large meteorites. On one of these they also measured 
the distribution of neon and argon. Contours of equal gas content indicate the 
approximate shape of the original meteorite before it entered the Earth’s atmo- 
sphere. By making use of a model for rare gas production from both primary and 
secondary radiation, estimates can also be made of the original size of the body and 
of the radiation dosage. 

Some very interesting results on the determination of cosmic ray exposure ages 
from the Helium-3 to Tritium ratio were described by J. Geiss (Bern). The 
Tritium content at the time of fall of a meteorite is a measure of the rate of produc- 
tion of Helium-3. (Some Helium-3 is produced directly by spallation, while the 
remainder is the daughter product of Tritium.) Consequently it is a simple matter 
to calculate an exposure age from the He*/T ratio, and there is no need to guess 
at the radiation intensity or self-shielding effects of the meteoritic material. One 
group of results contains ages between 20 and 24 million years, and the other group 
lies between oand 8 million years. Additional evidence from the radiogenic Helium-4 
age suggests that the latter group has lost Helium through being heated very late 
in its history. It is supposed that the 20~24 million year group is «vidence of the 
simultaneous formation of a number of meteorites by the breaking up of a much 


bigger body. 


5. Secular Variation and Palaeomagnetism 
5.1. Secular variation 


In his Chairman’s report Nagata presented isoporic charts for 1955-60, and 
pointed out the newly found intense foci in Antarctica, such large secular variation 
(200 y per year in Z) contrasting with that in the Arctic. Thellier described how, 
using palaeomagnetic data from archaeological specimens the secular variation in 
Europe has been traced back 2000 years. During this time the intensity of the 
Earth’s field has decreased by half. 


5.2. Origin of the geomagnetic field 


Most discussion centred on three papers by Rikitake, Lowes and Mathews, 
concerning theoretical representations of the magneto-hydrodynamics of the 
Earth’s core. Lowes and Mathews have independently investigated the behaviour 
of two coupled disk dynamos with the aid of computers, and have found that 
over a wide parametric range interaction between the dynamos produces sharp 
current reversals superimposed on a higher frequency oscillation. The reversals, 
of unpredictable frequency and duration, produce polarity reversals of the external 
magnetic field. Such models, though highly idealized, can present a coherent ex- 
planation of the geomagnetic field reversals indicated by palaeomagnetic data. 


5.3. Rock magnetism 


Emphasis was laid on two topics; chemical remanent magnetism (C.R.M.) and 
anisotropy of susceptibility. Professor Nagata began by drawing attention to the 
method of production of C.R.M. during a chemical transformation, the magnetic 
material changing from a “‘superparamagnetic” state into a ferromagnetic single 
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domain assemblage once the grains have grown large enough and then into 
multidomain assemblage when the grains grow too large to support single domains. 
Kume reported how the growth of C.R.M. in the goethite—hematite transforma- 
tion is affected by the application of pressure. ‘The new magnetism deviates 
from the direction of the applied field towards a direction perpendicular to that of 
the applied pressure. Howell noted in some rocks with chemical magnetization 
that the plane of maximum susceptibility coincides with the bedding plane which 
is also the plane of alignment of hematite crystals, and suggested that stress had 
played some part during their magnetization. Noltimier considered susceptibility 
anisotropy of red sediments as most likely due to crystallographic alignment. 

The use of alternating field demagnetization techniques was mentioned by a 
number of speakers. Mme Petrova demonstrated the relation between magnetic 
stability and the average coercive force of a rock, a general guide being that rocks 
with coercive force greater than 50 oersted are stable. 

On measuring remanence directions of dykes in the Canadian shield, Strangway 
noticed a correlation with their strike directions. He showed that if the internal 
demagnetizing field of a dyke is large enough it can deviate the remanence direc- 
tion from that of the ambient field. 


5.4. Palaeomagnetism 


New palaeomagnetic data were reported from many parts of the world, from 
Siberia to Antarctica. Amongst the most interesting were those for the Devonian 
from the western part of U.S.S.R. which are in substantial agreement with Euro- 
pean results. From a detailed study of a polarity reversal in a lava sequence in 
Kamchatka Sholpo was able to conclude that the Earth’s field had reversed its 
polarity and in so doing had decreased in magnitude as the vector rotated. Nagata 
reported the same phenomenon, indicating that the Earth’s field had decreased to 
a quarter, or less, of its normal magnitude as the reversal took place. Kropotkin 
demonstrated how the palaeoclimatic significance of the distribution of palaeozoic 
salt deposits in the U.S.S.R. was in agreement with the palaeomagnetic data, and 
Runcorn stressed in particular the need for a thorough testing of palaeomagnetic 
data in this way against palaeoclimatic evidence on a world-wide scale throughout 
the geological succession. Ina critical analysis of palaeomagnetic data Doell pointed 
out that although for the period since the Oligocene the results are consistently in 
favour of an axially dipolar geomagnetic field which every so often reverses its 
polarity, the results for older rocks can be plausibly interpreted in terms of a 
relatively rapid polar wandering. He maintained that the data are not yet suffi- 
ciently reliable for any case to be made either for or against the continental drift 
hypothesis, in contrast to a recent study by Blackett which favours continental 
drift. Doell voiced the general opinion in calling for more data, combined with a 
critical approach to the subject. 











Book Review 


Handbook of Geophysics 


Revised Edition. Multiple authors of Geophysics Directorate, 
United States Air Force. 


(84 in. X 11 in., 680 pp, Macmillan, New York, 1960, $15.) 


This is another of those “handbooks” which almost demand a trolley for their 
transport but, .his said, it is otherwise to be welcomed into the rapidly growing 
literature of geophysics. It originated in 1957 in the Geophysics Research Direc- 
torate as a handbook of geophysical data for internal consumption and for the U.S. 
Air Force’s engineers and designers and was not available to scientists generally. 
Having proved its worth it has been revised and enlarged and is now on sale to the 
public. In terms of dollars per pound weight it is not expensive but with better 
arrangement and lay-out, thinner paper, etc., it could probably have been reduced 
to about two-thirds its size and been the more handleable. 

The main point in reviewing a book of this kind is to describe the sort of 
information to be found in it and to assess its reliability. There are 21 chapters 
with titles (they should all be given if readers are to decide whether it is to be useful 
to them): 


Model Atmospheres Atmospheric Composition* The Ionosphere 
Temperature Atmospheric Electricity ‘Thermal Radiation 
Atmospheric Pressure Geomagnetism The Sun 
Atmospheric Density Meteors Cosmic Radiation 
Wind Terrestrial Surface Para- Jet Aircraft Condensa- 
Precipitation meters tion Trails 
Clouds Electromagnetic Wave Atmospheric Exploratory 
Propagation in the Devices 
Lower Atmosphere Acoustic Propagation of 
Visibility the Atmosphere. 


* Includes Aerosols. 


As these titles mainly imply, the handbook provides data primarily on the 
Earth’s atmosphere and on the skin of the Earth; there is nothing on the interior 
of the Earth or on the oceans except their surface properties. The chapter on ““Ter- 
restrial Surface Parameters” does, however, contain sections on, inter alia, Earth 
tremors—rather sketchy—and on geodesy, the latter being reproduced from part 
of W. A. Heiskanen’s Geodetic Data (McGraw-Hill, 1957). 

The first chapter, on Model Atmospheres, immediately raises the question of 
nomenclature which, with the increasing tempo of research on the upper atmo- 
sphere, should be so agreed as to avoid confusion among the many workers, with 
different disciplines, in the field. The Handbook lists most of those nomenclatures 
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in use until a few years ago but does not include the system based on thermal struc- 
ture which seems to be gaining the increasing assent both of conventional meteoro- 
logists and so-called aeronomers, namely: troposphere (no confusion), strato- 
sphere—from tropopause to temperature maximum (stratopause) at ca. 55 km, 
mesosphere—from stratospause to temperature minimum (mesopause) at 
ca. 85 km, and thermosphere—above mesopause. This is a simplification of the 
Chapman system and gets rid of the ill-defined level at ca. 25 km where Chapman 
placed the upper limit of the stratosphere and the lower limit of thy mesosphere. 
Under the more recent classification the stratosphere is the whole of the layer with 
a vertical temperature gradient markedly more stable than the troposphere below 
and the mesosphere above and is in line with the original meaning attached to 
“stratosphere” when sounding balloons were the only direct exploratory tool for 
high levels and rarely exceeded a height of z2okm or so. A lower and upper strato- 
sphere can be introduced into the system described above when a rough distinc- 
tion is to be made between the less and more stable parts respectively of the 
stratosphere. 

Apart from the oceanography and the Earth’s interior, referred to earlier, the 
data coverage is on the whole very comprehensive and up to date. Satellite and 
rocket data up to 1959 have for example been included. The engineering applica- 
tions of geophysical data are very much to the fore and this leads here and there to a 
welter of units but a scientist with preference for c.g.s. or c.g.s. derived units is 
not materially hindered. It also results in many tables and diagrams being included 
which will not be of interest to the geophysicist as such but, bulk apart, this is 
readily to be tolerated for the wealth of material which the geophysicist will be glad 
to have in so accessible aform. In most sections it is made clear that data gathered 
from one part of the world (very often the U.S.) may not be representative of other 
parts and that where knowledge is sparse, e.g. in the higher atmospheric levels for 
many variables, only a first approximation is provided and that often by little 
more than inspired guesswork. There is a seeming reluctance to show or use any 
theoretical relationships which can be so important in welding geophysical data into 
a self-consistent whole, though such relationships must have formed the basis of 
much that appears here. This again no doubt derives from the volume being 
addressed to engineers as much as to scientists but the latter are thereby denied the 
means of assessing reliability in many cases. A simple example arises in relation 
to the first four chapters in which the hydrostatic equation puts a remarkable straight- 
jacket (with the equation of state) on the variations of temperature, pressure and 
density. ‘The reviewer has found no mention of this equation. There is, however, 
a very substantial and, where tested by this reader, seemingly comprehensive set 
of references to the sources of data included in the volume. 

An interesting hangover (I presume) from the first restricted edition is that here 
and there where only a first approximation has been provided for the variation of 
some property in space and time the reader is invited to request additional informa- 
tion or even special studies for his own needs. What a handsome service G.R.D. is 
offering us, or are they? 

There are of course many items of knowledge which a geophysicist might hope 
to find within these covers but will not. The reaction of most, however, will more 
likely be one of amazement and gratitude that so much patient effort has been 
expended in amassing so much information and presenting it so clearly—many of 
the diagrams are indeed given on an unnecessarily large scale. Of the reviewer’s 
own noted omissions one worth mentioning here is the absence, in the otherwise 
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excellent chapter on Electromagnetic Wave Propagation in the Lower Atmosphere, 
of the effects of refraction on optical transmission—a matter surely of wide concern. 
To offset this mild ‘criticism the reviewer would, on a very personal basis, like to 
express a great admiration for the chapter on Thermal Radiation and, on behalf of 
all experimental geophysicists, welcome particularly the penultimate chapter on 
instrumental techniques and performance. Other users will no doubt select other 
chapters as particularly valuable or well-conceived and, perhaps, find others less 
than adequate. No one man today is able to assess the balance and validity of 
material over the whole range covered by this volume, but it is at least evident that 
a valuable service has been provided the geophysicist and one can look forward to 
later editions in which imperfections and omissions, as they appear, will be recti- 
fied. There is an 11-page index and the volume is produced to withstand much 
handling. This is a “must” for every geophysical library and for the personal 
shelves of all who can afford it. 


P. A. SHEPPARD 





On Mr King-Hele’s Theory of the Effect of the Earth’s 
Oblateness on the Orbit of a Close Satellite* 


P. J. Message 


(Received 1960 February 8) 


Summary 


The relations are obtained between « and Q, the parameters 
defining King-Hele’s plane 7, and the osculating inclination and 
longitude of node. The secular motion of the node is thereby obtained 
in terms of mean osculating elements in a form comparable to that arising 
in Brouwer’s treatment of the canonical equations. A modification is 
made to King-Hele’s theory which removes the singularities in the 
Je terms of dQ/dis, and it is also found that there are no terms of order 
Je? in « or Q. 


* The full text of this paper is to be found in Mon. Not. R. Astr. Soc., 121, 1 1960. 
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Geophysical Discussions of the Royal Astronomical 
Society 


The following Geophysical Discussions have been arranged for 1961 and will 
take place in the Society’s rooms at Burlington House, Piccadilly, London W.1 
at 16" 15™ (tea at 15" 45™): 


1961 January 27.— The geological time scale - 
1961 February 24.—Upper atmospheric ionization and aurorae 
1961 March 24.—Geodetic uses of artificial satellites. 


No meeting for the discussion of papers has been arranged for this session. 
Fellows are reminded that they may ask to read papers at any ordinary meeting of 


the Society. 
A. H. Cook 


Geophysical Secretary 





ERRATUM 


Geophys., ]., 3, No. 1, 1960: 
G. Bomford, The figure of the Earth—its departure from an exact spheroid, 
Pp. 92, line 16: 
For (i) +0-008/4/(L) metres in Alpine country, 
read (i) +0°008y/(/L) metres in Alpine country. 
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| NOTICE 

“The Council has pleasure in announcing publication, early in 1961, of an 
extra volume of the Geophysical Jowr=c! (Volume 4), in honour of the seventieth 
birthday of 


Sir Harold Jeffreys 
The Geophysical Editors have obtained the co-operation of many of 
Sir Harold’s distinguished colleagues and pupils and believe that the contents of 
this special volume will represent the most advanced work in a number of fields 
in which Sir Hzzold has long been interested. 


Contrisutors INCLUDE: 
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Bi Sayad M. Hassan 
K. Hidaka 

U. Hochstrasser 

J. E. Jackson 

J. A. Jacobs 

H. Jarosch 

N. Jobert 

D. G. King-Hele 

L, Knopoff 

R. L. Kovach 

¥. R. Lapwood P. L. Willmore 
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Z. 
M. 
F. 
E.C 
K. E. Bullen 
P. 
A. 
A. 
F. 
FF. 
A. 


This volume, suitably bound in cloth boards, will also be available as a special 
publication. 


Price £4 (US$ 12) 


Orders, with remittance, should be sent to: 
Assistant Secretary, Royal Astronomical Society, Burlington House, London, W.1. 





SPECIAL NOTICE TO FELLOWS 


” Fellows of the Society may obtain this bound volume of the Geophysical Journal 

; ially reduced price of {1 108. od. (US$ 4.50). It is particularly requested 

that ‘ellows wishing to receive a copy should send their orders, with remittance, 
as soon as possible to: 


- Assistant Secretary, Royal Astronomical Society, Burlington House, London, W.1. 
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